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ABSTRACT 

We report extensive radial-velocity measurements of the two giant components of the detached, 104- 
day period binary system of Capella. Our highly accurate three-dimensional orbital solution based 
on all existing spectroscopic and astrometric observations including our own yields much improved 
masses for the primary and secondary of 2.466 ± 0.018 Mq and 2.443 ± 0.013 Mq, with relative errors 
of only 0.7% and 0.5%, respectively. The mass ratio is considerably closer to unity than previously 
believed, which has an impact on assessing the evolutionary state of the system. Improved values are 
presented also for the radii (11.87±0.56 Rq and 8.75±0.32 Rq), effective temperatures (4920±70 K 
and 5680 ± 70 K), and luminosities (79.5 ± 4.8 Lq and 72.1 ± 3.6 Lq). The distance is determined to 
be 13.042 ± 0.028 pc, based on the accurate orbital parallax. The projected rotational velocities and 
individual rotation periods are also known. Capella is unique among evolved stars in that, in addition 
to all of the above, the chemical composition is known as well. This includes the overall metallicity 
[m/H], the carbon isotope ratio ^^G/^^G for the primary, and the lithium abundance and C/N ratios for 
both components. We present new or revised values for some of these. The latter three quantities are 
sensitive diagnostics of evolution, and change drastically for giants as a result of the deepening of the 
convective envelope during the first dredge-up. The secondary is crossing the Hertzprung gap, while 
the primary is believed to be in the longer-lived core-helium burning phase. Previous studies using 
only the masses, temperatures, and luminosities have found good agreement with stellar evolution 
models placing the primary in the clump. Here we compare all of the constraints simultaneously 
against three sets of current models. We find that they are unable to match all of the observations 
for both components at the same time, and at a single age, for any evolutionary state of the primary. 
This shows the great importance of chemical information for assessing the evolutionary state of giant 
stars. A comparison with models of tidal evolution yields similarly disappointing results, when tested 
against the fact that the orbit is circular, the primary is rotating synchronously, the secondary ~12 
times faster than synchronous, and the spin axes are apparently aligned with the axis of the orbit. 
When confronted in detail, our understanding of the advanced stages of stellar evolution is thus still 
very incomplete. 

Subject headings: binaries: general — binaries: spectroscopic — stars: abundances — stars: evolution 
— stars: fundamental parameters — stars: individual (Capella) 



1. INTRODUCTION 

Our knowledge of stellar evolution relies heavily on 
constraints provided by measurements of the physical 
properties of stars, particularly the mass. Eclipsing bina- 
ries have been our main source of accurate information 
for stars on the main sequence (masses, radii, effective 
temperatures, luminosities). The more advanced stages 
of evolution, after the central hydrogen is exhausted, 
are much less well understood than the main-sequence 
phase primarily because of poorer or less stringent con- 
straints. Accurate masses are rarely available for normal 
giant stars, for practical reasons that are well known: 
it requires them to be in binaries, and it requires the 
orbits to be wide enough so that the components are de- 
tached (i.e., that they do not interact with each other in 
the sense of mass transfer). This in turn calls for sus- 
tained observations over long periods of time, often sev- 
eral years. The radial-velocity amplitudes are typically 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden St., 
Cambridge, MA 02138, USA; e-mail: gtorres@cfa.harvard.edu 

^ Instituto de AstroHsica de Andalucfa, CSIC, Apartado 3004, 
18080 Granada, Spain; e-mail: claret@iaa.es 

^ School of Earth and Space Exploration, Arizona State 
University, P.O. Box 871404, Tempe, AZ 85287; e-mail: 
patrick.young.l@asu.edu 



low, making the masses more difficult to measure accu- 
rately. Eclipses are also much less likely, robbing us in 
most cases of the possibility of determining the absolute 
radii directly. 

While precious few normal giants or subgiants 
are known to be in eclipsing binaries (examples in- 
clude AI Phe, TZ For, and the system OGLE 
0510 1 9.64-685812.3 in the LM C; I Andersen etldlfigSSL 
119911: iPietrzvhski et all l2009t l. non-eclipsing systems 
amenable to both spectroscopic and astrometric stud- 
ies enabling accurate dynamical mass determinations are 
somewhat more common. A very prominent example of 
this class of objects is Capella [a Aurigae, HD 34029, 
HR 1708, HIP 24608, G8 III-I-GO III, P = 104 days), the 
6th brightest star in the sky {V = 0.07), and the first bi- 
nary for which an astrometric orbit was determined inter- 
ferometrically (Anderson 1920; McrriU 1922). Although 
other astrometric measurements have been gathered for 
Capella over the decades using a variety of techniques, 
those pioneering observations with the original Michelson 
interferometer on Mount Wilson have remained a criti- 
cal component of the astrometric solution of the orbit, 
which until about 15 years ago contributed the dominant 
shar e of the uncerta i nty in the resulting masses. At that 
time iHummel et~all () 19941) used the Mark HI interfer- 
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ometer, also on Mount Wilson, to obtain measurements 
that are an order of magnitude more precise. As a re- 
sult, our knowledge of the m asses is now l i mited by the 
best-available spectroscopy () Barlow et al.l |1993[ ) rather 
than the astrometry, a somewhat embarrassing situation 
given the brightness of the object and its 100-year obser- 
vational history. The masses are currently estimated to 
be 2.69 ±0.06 M© and 2.56 ±0.04 Mq for the cooler pri- 
mary and the secondary, respectively. Although these are 
seemingly quite precise (~2% relative errors), the accu- 
racy is just as important, especially for giants. The pre- 
dicted properties of stars from stellar evolution models 
in these rapid evolutionary stages are extremely sensitive 
to mass, and a 2% error makes a much larger difference 
than on the main sequence. 

Systematic errors in the radial velocity measurements 
for the rapidly-rotating secondary of Capella are not easy 
to avoid, and as a result there have been persistent con- 
cerns in the literature about the accuracy of the mass 
ratio, which is very close to unity. This quantity has a 
significant impact on the precise evolutionary state in- 
ferred for the components. One of our motivations here 
is thus to provide a new, high-quality set of spectroscopic 
observations that more nearly matche s the p recision of 
the astrometric data of iHummel et al.l (|1994f). We focus 
especially on the accuracy of the mass determinations, 
for the purpose of comparing with state-of-the-art stel- 
lar evolution models. Considerable effort is therefore in- 
vested in the inter-comparison of all available astrometric 
and spectroscopic observations, in order to properly un- 
derstand the systematics. An additional goal of this work 
is to carry out a comprehensive analysis of all the data 
relevant for the determination of the evolutionary state of 
Capella, which has been a subject of debate for decades. 
The secondary is clearly crossing the Hertzprung gap, 
but the primary has been suggested to be cither in the 
helium-burning clump or on the first ascent of the giant 
branch, usually based on partial information. This sys- 
tem is perhaps unique in that, in addition to the masses, 
effective temperatures, and luminosities that have been 
used previously for that purpose, a wealth of other infor- 
mation is available. This includes direct measurements 
of the angular diameters, various activity indicators in 
the optical, ultraviolet, and X rays, the projected ro- 
tational velocities as well as the rotation periods, the 
overall metallicity, and particularly the surface lithium 
abundance of both stars, the ^^C/^^C isotope ratio for 
the primary star, and the C/N abundance ratios. Chem- 
ical indicators such as these are crucial diagnostics of 
evolution because they change significantly in the giant 
phase, mainly during the first dredge-up. We bring all of 
these to bear here, for the first time. We wish to examine 
the effects of convection prescriptions in the models (mix- 
ing length, overshooting) as well as rotation, which has 
not previously been considered for this system. Further- 
more, Capella is an important point of comparison with 
tidal evolution theory for evolved stars. This is because 
the primary has its rotation synchronized with the or- 
bital motion while the secondary rotates 12 times faster 
than synchronous, despite the nearly identical masses, 
which arc within 1% of each other. Thus another goal 
of this work is to use our newly derived accurate dimen- 
sions for the components to carry out a detailed check 
of current tidal theories and gauge our understanding of 



these processes. 

Our new spectroscopic observations are presented in 
§[5J along with all historical radial-velocity measure- 
ments. The astrometric observations are discussed in §|3l 
and include long-baseline interferometry, speckle interfer- 
ometry, and Hipparcos measurements. Our simultaneous 
orbital solution for Capella based on these two kinds of 
data is documented in §|4l with particular care given to 
possible systematic errors that might bias the masses. In 
§[5] and §|n]we collect all the information available on the 
relative brightness of the components and the angular 
diameters, which we use later to estimate effective tem- 
peratures and absolute radii. The chemical composition 
of Capella is discussed in §[3 and proves to be of criti- 
cal importance for the analysis. The physical properties 
of the stars are then described in §|8l We present a de- 
tailed comparison of the absolute dimensions of Capella 
with stellar evolution theory in §|9l focusing on the de- 
termination of the evolutionary state of the stars. In the 
same section we carry out tests of current tidal theories. 
Finally, §[TU] summarizes our main conclusions. Two ap- 
pendices collect notes of interest on the astrometric ob- 
servations for the benefit of future users, and a third con- 
tains a discussion of the coronal abundances of Capella 
that support the photospheric determinations. 

2. SPECTROSCOPIC OBSERVATIONS 

The rich history of the radial- velocity measurements of 
Capella began more than a century ago with the pioneer- 
ing efforts by the Gree nwich and Potsdam astronomers 
(|GiIilll891l: rVogellll891h . and has b een recounted prev i- 
ously by other authors (see, e.g., iBarlow et al.l Il993l ). 
More than a dozen semi-independent spectroscopic or- 
bital solutions have been reported over the decades, 
based on data sets of greatly varying quality. Despite 
the brightness of the object, the radial velocity measure- 
ments of the rotationally broadened secondary compo- 
nent are not particularly easy, and have not always been 
possible in the past. As a result, there has been consider- 
able debate about the mass ratio (see, e.g.. lBatten et al.l 
|1991|) . which our analysis conclusively shows is near 
unity. Nevertheless, some of these historical data are 
still of value to improve the orbital period, so we de- 
scribe them in some detail below and make use of them 
later. Our main observational contribution here is a large 
new set of high-quality velocity measurements for both 
components that provides more than a two-fold improve- 
ment in the precision of the velocity semi-amplitude of 
the primary, and a four-fold improvement for the sec- 
ondary comp a red to the best existing determinations, by 
IBarlow et all (|1993f ). These new data are described first. 

2.1. New radial-velocity measurements 

Spectroscopic observations of Capella were conducted 
at the Harvard- Smithsonian Center for Astrophysics 
(CfA) using the 1.5m Wyeth refiector at the Oak 
Ridge Observatory (Harvard, Massachusetts), beginning 
in 1996 October and continuing through 1999 Novem- 
ber. An echelle spectrograph with a photon-counting in- 
tensified Reticon detector (Digital Speedometer: Lathainl 
[T9851IT991 was used to record a single 45 A echelle or- 
der centered at a wavelength of 5188.5 A, featuring the 
gravity-sensitive lines of the Mg I b triplet. The re- 
solving power provided by this setup is A/AA w 35,000. 
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One additional observation was gathered with a nearly 
identical system on the 1.5m Tillinghast reflector at the 
F. L. Whipple Observatory (Mount Hopkins, Arizona). 
Nominal signal-to-noise ratios per resolution element of 
8.5 kms~^ range from about 20 to 90, although for values 
much higher than 50 the limit is set by systematics from 
flat-fielding and not photon noise. With the inclusion 
of two archival observations made in 1986 February and 
March with the instrument at Oak Ridge, the total num- 
ber of usable spectra is 162, collected over an interval of 
13.8 years. 

Radial velocities for both stars were derived using 
TODCOR, a t wo-dimensional cro s s-corr elation technique 
introduced by IZucker fc Mazehl (|1994D . This method 
uses two templates, one for each component of the bi- 
nary, which we selected from a large library of synthetic 
spe ctra based on mode l atmospheres by R. L. Kurucz 
(see lLatham et al.|l2002D . These templates have been cal- 
culated for a wide range of effective temperatures (Toff), 
surface gravities (log g) , rotational velocities {v sin i when 
seen in projection) , and metallicities ([m/H]). Following 
[Torres et al.l ()2002[ ) the optimum templates for Capella 
were determined by means of extensive grids of cross- 
correlations with TODCOR, seeking to maximize the av- 
erage correlation weighted by the strength of each expo- 
sure. The surface gravities were held fixed at preliminary 
values of logg — 2.5 and 3.0 for the primary and sec- 
ondary, and the metallicity was initially assumed to be 
solar. As a result of this optimization we obtained by in- 
terpolation effective temperatures of — 4900 ± 100 K 
and = 5710 ± 100 K for the primary and secondary, 
respectively, along with projected rotational velocities of 
^;Asinj = 6.5 ± 1.0 kms~^ and ?;Bsinj — 36.0 ± 1.5 
kms~^. Strictly speaking, the latter values are a mea- 
sure of the total broadening of the spectral lines. We 
discuss these measures and compare them with others in 
§[8l We repeated these determinations assuming metal- 
licities of [m/H] = —1.0, —0.5, and -1-0.5, but we found 
the average correlation values to be slightly lower than 
with solar metallicity, indicating a poorer match to the 
observed spectra. The templates adopted here are the 
ones in our library with parameters nearest to the val- 
ues above: Teff = 5000 K and 5750 K for the primary 
and secondary, and rotational velocities of 6 kms~^ and 
35 kms""'^, respectively. The radial velocities derived in 
this way have internal errors averaging 0.5 km and 1.0 
kms~^ for the primary and secondary, but vary individ- 
ually depending on the signal-to-noise ratio. In addition 
to the radial velocities, our spectra yield the light ratio at 
the mean wavelength of our spectra, ^b/^a = 1.48±0.05. 
The difference in line blocking between the components 
has been accounted for, so that this represents a true flux 
ratio rather than a ratio between the continuum levels. 
The hotter star is thus brighter at ^5200 A. 

The stability of the zero-point of the CfA velocity sys- 
tem was monitored by means of exposures of the dusk 
and dawn sky, and small systematic run-to-ru n correc- 
tions were applied in the manner described by iLathanJ 
(|1992| ). The zero point of the native CfA velocity sys- 
tem based on synthetic templates is very close to the 
absolute frame as defined by extensive observations of 
the minor planets in the solar system. The correc- 
tion required to place our radial velocities on this abso- 
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Fig. 1. — Systematic errors in the raw TODCOR velocities for 
Capella as a function of orbital phase and velocity. Filled circles 
correspond to the primary, open circles to the secondary. These 
differences have been applied to the raw velocities as corrections 
(see text). 



lute frame is -1-0.139 km s~ ^ (|Stefanik. Latham fc Torrei 
119991 : [Latham et al.ll200^ . and has not been applied to 
the measurements listed below. 

One of the main advantages of TODCOR compared 
to conventional one-dimensional cross-correlation tech- 
niques is that it greatly reduces the systematic errors in 
the radial velocities caused by line blending that have 
affected many of the previous studies of Capella (see 
S I2.2p . Nevertheless, because we are concerned in this 
work with the accuracy of the velocities as much as their 
precision, we investigated possible systematic effects that 
may result from residual blending in our narrow spectral 
window, and from shifts of the spectral lines in and out 
of this window as a function of orbital phase. Previous 
experience with similar material has shown that these ef- 
fects are sometimes significant, and must be examined on 
a case-by-case basis and corrected if necessary (see, e.g., 
[Torres et al.lll997l I200C ) . We perforni e d num erical simu- 
lations as described bv lLatham et aTl (jl996l ) to evaluate 
these effects. Briefly, we generated synthetic composite 
spectra matching our observations by combining the pri- 
mary and secondary templates used above (including ro- 
tational broadening), shifted to the appropriate velocities 
for each of the exposures as predicted by a preliminary 
orbital solution, and scaled to the observed light ratio. 
These synthetic observations were then processed with 
TODCOR in exactly the same way as the real spectra, 
and the resulting velocities were compared with the input 
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shifts. The differences are shown graphically in Figure [U 
as a function of both velocity and orbital phase. The 
systematic pattern is obvious, and the individual differ- 
ences can reach ±1 kms~^, which is relatively small in 
absolute terms but significant compared to the internal 
errors. We therefore applied these differences as correc- 
tions to the raw velocities. The effect on the primary 
semi-amplitude Kp^ is negligible, but the change in 
is ~0.5%, which translates into a non-negligible change 
in the derived masses of about 1% for the primary and 
0.6% for the secondary. The final velocities in the he- 
liocentric frame are given in Table [l] and include these 
corrections. Similar adjustments based on the same sim- 
ulations were applied to the light ratio, and are already 
included in the value reported above. 

Preliminary single-lined orbital solutions performed 
separately on the primary and secondary velocities indi- 
cated a slight difference in the center-of-mass velocities, 
7, of about 0.27 ± 0.08 kms^^, with the secondary value 
being lower. Primary/secondary velocity differences con- 
siderably larger than this are not uncommon in studies of 
double-lined eclipsing binaries. This difference Aab per- 
sisted in our global solution described later. Although 
it is very small in absolute terms (only about half of the 
typical uncertainty in our primary velocities), it is statis- 
tically significant due to the large number of observations 
in the fit. Because it may affect the absolute masses of 
Capella at some level, we have explored possible reasons 
for this shift. One is the differential gravitational redshift 
between the stars, given that our synthetic templates do 
not account for this. Estimates based on preliminary 
values of the masses and radii of the components indi- 
cate that this effect is 0.046 kms~^, but it goes in the 
wrong direction to explain Aab, i-e., the secondary red- 
shift is larger. It is also possible there are shift s due 
to large- scale convective motions ([Schwarzschildl 119751 : 
iPorter fc Woodwardll2000[ l that could be different in the 
two stars, but these are not well characterized for giants. 
Given that the stars are observed to be active, another 
possibility is the presence of spots on one or both com- 
ponents (particularly on the rapidly rotating secondary) , 
which can affect the velocities. A perturbation of this 
nature was in fact pointed out bv lHummel et all ()1994l ) 
for their interferometric visibilities of Capella (see §|4]). 
Unfortunately our time sampling is inadequate to study 
this in more detail, but unless the spots are very long- 
lived we would expect the effect to average out to some 
extent over the interval of our observations. A fourth 
possibility that cannot b e completely rule o ut is tem- 
plate mismatch (see, e.g.. iGriffin et al.ll2000[ ). We have 
made every effort here to use templates that maximize 
the average correlation for all our spectra, and small dif- 
ferences with the true values of Tcff, tisinz, logg (which 
we estimate below to be \ogg = 2.68 for the primary and 
2.94 for the secondary) or metallicity compared to what 
we have assumed should not have a significant effect on 
the velocities, in our experience. However, line broaden- 
ing from micro- or macro-turbulence in Capella could be 
somewhat different from what is assumed in our library 
of synthetic spectra (microturbulence = 2 kms~^ and 
macroturbulence Crt = 1-5 kms~^), although this is 
unlikely to affect the secondary much due to the over- 
whelming effect of rotational broadening in that star (36 
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kms ). We discuss this effect further in § 19.21 in con 
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Fig. 2. — CfA radial-velocity measurements for Capella (filled 
circles for the primary, open for the secondary) along with the 
curves computed from our combined orbital solution described in 
§|4] The dotted line represents the velocity of the center of mass. 
O — C residuals for the primary and secondary are shown in the 
bottom panels. 

nection with the accuracy of the vsini measurements. 
In the absence of a definitive explanation, we have cho- 
sen below to correct for the primary/secondary shift by 
solving for the offset and applying it to the secondary 
velocities. Not correcting for the shift would affect the 
semi-amplitudes at the level of 0.03% for the primary and 
0.14% for the secondary, and the final masses at the level 
of 0.31% and 0.22%, which correspond to less than half 
of the formal uncertainties in our final determination of 
those quantities (see §[8|). 

Figure [2] displays the CfA observations along with an 
orbital solution described later, as well as the O—C resid- 
uals. Those of the secondary show a small residual pat- 
tern which we believe to be of a similar nature as the 
Aab shift discussed above. We return to this later in 
connection with the orbital solution. 

2.2. Historical radial- velocity measurements 

The discovery of the binary character of Ca pella 
was aiinounce d independently bv ICampbelll ()1899f ) and 
iNewalll ()1899[ ) from photographic spectra taken at Lick 
Observatory and Cambridge Observatory (England) , re- 
spectively. Both investigators noted the composite na- 
ture of the spectrum, but published velocities only for the 
component Newall referred to as being of "solar type" . In 
our nomenclature this is the cooler, slightly more mas- 
sive star we refer to as the "primary" (star A), which 
has relatively sharp lines. The other star (of "Procyon 
type", "secondar y " , or star B) has much more diffuse 
lines. I Campbell (|1901f ) reported only that the veloc- 
ity of the secondary varies between —3 kms^^ and -t-63 
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kms^^. His 31 measurements for the primary are of ex- 
cellent quality ((TRv ~ 0.8 kms~^), and were used by 
iReesd (|190Clf ) to establish t he first reliable spectroscopic 
orbit. The measurements bv lNewaUl ()1900( ) are somewhat 
poorer (ctrv ^2.1 kms^^), but still potentially useful. 

The first published measure ments of the secondary 
velocity appear to be those by iGoosI (I1908D . who suc- 



ceeded in detecting it in 19 of his 35 photographic plates 
taken with a 0.3m refractor at Bonn. All plates yielded 
good measurements for the primary. Further veloc- 
ities for both stars we re repo r ted b y ISanfordI ({19221 ) 
from Mount Wi l son, a nd lStruvd ()1939l ) from Babelsberg. 
IStruve fc Kilbyl (|1953l ) published a further series of veloc- 
ities from Mount Wilson and Lick, though only for the 
primary star. Measurements of the velocity difference 
between the components made at the Dominio n Astro- 
physic al Observatory (DAO) were reported by IWrightj 
(|1954 ). along with a careful study of the secondary spec- 
trum and the relative brightness of that star compared 
to the primary. The velocity measurements are averages 
from plates taken at similar orbital phases, so unique 
dates cannot be assigned to them and for this reason 
we do not use these data here. The brightness ra tio as 
well as the mass ratio estimate by IWrightj (|1954|). de- 
rived by adopting the primary orbit from* Struve fc Kilbvl 
(jl953|), were quite influential over the following decades, 
although the light ratio is now known to be incorrect 
(or at least misl e ading ; sec §[5] and §|S]). More recently 
iBatten fc Breed (|1975l ) published 18 velocities for the 
primary from plates obtained at DAO, rnost of which 
were later re- measured by Batt en et al.l (|1991f ). super- 
seding the original determinations. Further measure- 
ments for both compone nts obtained at th e Fick Obser- 
vatory were reported by IShen et aD (jl985f ) . These data 
were in turn sup erseded and significantly expanded by 
iBeavers fc Eittej |i986), who published the largest set 
of velocities for Capella aside from our own. Additional 
measurements of t he pr imary only were obtained by 
IShcherbakov et ahl (jl990D . including a set based on pho- 
tospheric lines and another set from the chromospheric 
He I A10830 line. We do not use the latter because they 
may not correspond to the true center of mass of the 
star, nor do we con sider a similar list of velociti es for 
both components bv iKatsova fc Scherbakovl ()1998l ). also 
from the He I A10830 line. Finally, high-quality measure- 
ments for both stars from the M cDonald and Kitt P eak 
Observatories were published bv lBarlow et al.l ()1993f ). 

The sources above represent the most important veloc- 
ity data sets published for Capella in the century since 
its discovery as a binary. Even though some of them 
may have considerably larger uncertainties (less weight) 
than the CfA velocities, in principle there is no rea- 
son why they cannot be properly combined with ours to 
strengthen the solution, which is our goal in §|4l A num- 
ber of smaller lists of less than half a dozen measurements 
each have also appeared over the decades, but are ignored 
here for being much less significant and more difficult to 
use because of the poorly determined zero-point offsets. 
The richer sources are summarized in Tabled where the 
last entry corresponds to our own contribution. 

The potential usefulness of these historical data sets 
depends on whether they can be shown to be sufficiently 
free from systematic errors. To this end, we have exam- 
ined each of the sources by computing separate orbital 
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Fig. 3. — Radial-velocity semi-amplitudes for the primary and 
secondary of Capella (filled and open circles, respectively) through- 
out the last century, based on our solutions in Table [2] 



solutions from the original velocities with the same fit- 
ting code, and comparing them to one based on the CfA 
data. These solutions can be found in Table [2] We list 
also the number of observations, their time span, and 
the root-mean-square (RMS) residual from the fit in each 
case, which is representative of the typical error of the 
velocities. Because of the limited duration of some of 
these studies, the period has been held fixed at the value 
P — 104.022 days determined from a preliminary fit to 
our own observations, and the orbit has been assumed to 
be circular. The center-of-mass velocities 7 in the sec- 
ond column show that there are occasional differences 
in the instrumental zero points, although these can eas- 
ily be corrected in a combined solution by solving for 
the offsets simultaneously with the other adjustable pa- 
rameters. The same holds for the primary/secondary 
offsets listed in the third column (see also ii l2.1|) . which 
have been set to zero when a preliminary fit indicated 
the shift was not significant. We are more concerned 
here with the velocity semi-amplitudes Ka and Kb, 
which determine the masses of the components. Exclud- 
ing the two data sets that rely on chromospheric line s 
(IShcherbakov et al.lfT990t IKatsova fc Scherbakovlll998f l. 
the primary semi-amplitudes from all the others agree 
reasonably well with ours. The only exception is the 
data set bv I Sanford) (fl922l ). which is also the smallest. 
The secondary semi-amplitudes, however, show signifi- 
cant systematic differences with the CfA value of Kb, 
and if we restrict ourselves to velocities based on photo- 
spheric lines, there appears to be a trend of decreasing 
amplitudes as a function of time over the last century, 
leading up to our own determination (see Figure [3]). We 
suspect these differences have to do with systematic ef- 
fects associated with line blending and the difficulty of 
measuring the broad spectral features of the secondary, 
particularly in the older studies, a problem that has been 
pointed out repeatedly over the years. For this reason, 
we have chosen not to use any of the historical secondary 
velocities here, relying only on our own. The primary ve- 
locities, on the other hand, appear reason ably free from 
systematics (save those of ISanfordI I1922L which we ex- 
clude), and add up to more than twice the number of 
our own observations although the combined weight is 
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Phase 

Fig. 4. — Historical radial- velocity measurements for the primary 
component of Capella, along with the curve computed from our 
combined orbital solution described in §|4] Individual data sets 
have been placed on the same zero-point as the CfA velocities by 
applying the offsets described in §|3] The dotted line represents 
the velocity of the center of mass. O—C residuals are shown in the 
bottom panel. 

actually ~50% lower. 

All of these velocities are listed in Tablc[3]on their orig- 
inal scales (i.e., without the application of any offsets, to 
be described below).'* Individual uncertainties are de- 
scribed in §131 These observations are shown graphically 
in Figure |4l along with the same curve for the primary 
from Figure [21 

3. ASTROMETRIC OBSERVATIONS 

Soon after the discovery that Capella is a spectroscopic 
binary, some of the most skilled double-star observers of 
the day attempted to resolve the pair at the times pre- 
dicted from the spectroscopic orbit t o be the mo s t favor- 
able, but they wer e unsuccessful ((HussevI Il900t lAitkenI 
[ToMlHusse-ilflQnl . At about the same time, an intrigu- 
ing series of visual measurements was made at the Green- 
wich Observatory that did appear to barely resolve the 
object: the observers reported elongated images with the 
28-inch refractor. Systematic measurements of the posi- 
tion angle of the binary from these elongated images were 
carried out over an interval of about two years, and rough 
estimates of the angular separation were also made on a 
few occasions even though it was far smaller than the 
resolving power of the telescope. However, these obser- 
vations were never confirmed and have been called into 
question, so we do not use them in our analysis. Never- 
theless, a number of fascinating aspects of this puzzling 
data set are worth noting and are described in more de- 
tail in Appendix [Al 

It was not until 1919 that Capella was resolved in 
earnest, with the 6-m baseline Michelson interferometer 

4 The velocities by ICampbelll ||1901|) used here include small 
adjustments later determined bv lCampbell fc Moord (| 19281 ) to be 
required in order to place them on the scale of the homogenized 
catalog of 1896-1921 Lick velocities they published. These adjust- 
ments are specific to the person who measured the plates, and for 
this case are +0.1 kms~^ (Campbell) and —0.3 kms~^ (Wright). 



on the 100-inch te lescope on Mount Wilson (jAndersonl 
[19201: lMerrilil[l92l . These pioneering observations are 
of high quality and internal consistency and have since 
been used in nearly all of the astrometric orbital solu- 
tions published for the system. They are valuable be- 
cause of the extended time baseline they afford. We 
incorporate them into our own analysis as well, al- 
though they do contain some systematic errors that 
we addres s later. Except for two more recorded at- 
tempts by IWilsonI (|1939l . Il941f) to resolve the pair vi- 
sually, almost 50 years elapsed until the next astromet- 
ric observation s were made at Pulkovo Observatory by 
iKulaginl ()1970l ). with a similar interferometer also us- 
ing a 6-m baseline. Additional long-baseli ne interfero- 
metric observations have b een reported b y Blazit et al.l 
(fr977a) fbas ehne 12-20 mV.lKoechlin et al.' (1979) (base- 
line 13.8 m). iKoechlin et al.l (|1983D fbaseline 5.5-35 m), 
[Baldwin et al.l (|1996f ) (three-element Cambridge Opti- 
cal Aperture Synthesis Telescope, COAST, using base- 
lines up to 6.1 m), and more recently bv iKraus et al.l 
((2005 ) (three-element Infrared Optical Telescope Array, 
IOTA, using baselines up to 38 m). By far the most 
pr ecise interferornetric o bservations of Capella are those 
of [Hummel et al.l (|1994[ ) with the Mark III interferome- 
ter on Mount Wilson, using baselines of 3.0 to 23.6 m. 
These observations improved the uncertainties in both 
the position angle and the angular separation by about 
an order of magnitude compared to previous measures. 
They are also the only ones, aside from those obtained in 
1919-1921, that provide full phase coverage of the orbit. 
All interferometric measurements are listed in Table 2] 
(for those published in polar coordinates) and Table \5\ 
(measures published in Cartesian coordinates). 

Because of its brightness and convenient angular sep- 
aration, Capella has served for decades as an ideal 
calibration object for long-baseline interferometry, and 
has been referred to as "an interferometrist's friend" 
([Hartkopf et al.ll200lD . Be ginning in the 1970s Capella 
was observed also with the speckle interferometry tech- 
nique by a large number of investigators. Though typ- 
ically less precise that the long-baseline interferometry 
results, these measures are still useful and are folded into 
our solution below. They are collected in Table [6l 

Capella was a.lso a target of the Hipparcos mission 
([Perrvman et al.|[i997l ). It was observed under the desig- 
nation HIP 24608 a total of 43 times over a 3-yr interval 
(1990.08-1993.15), corresponding to nearly 11 orbital cy- 
cles of the binary. Each measurement consisted of a one- 
dimensional position ('abscissa', v) along a great circle 
representing the scanning direction of the satellite, tied 
to an absolute frame of reference known as the Inter- 
national Celestial Reference System (ICRS). The typical 
precision of these measurements is about 2.3 milli-arc 
seconds (mas) for Capella. The data were used by the 
Hipparcos team to solve for the five basic astrometric pa- 
rameters of the star, which are the position and proper 
motion components, and the parallax. Although the 
satellite measurements did not actually resolve the pair 
(separation ~56 mas), the motion of the center of light is 
large enough that it was clearly detected. Consequently, 
extra terms were added during the original reductions 
by the Hipparcos team to model this orbital motion and 
avoid biases. Several of th e orbital elements we re held 
fixed at the values from the iHummel et al.l (|1994l ) study 
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(period, epoch of nodal passage, inclination angle, posi- 
tion angle of the node), and the orbit was assumed to be 
circular. The semimajor axis of the photocentric motion 
reported in the catalog is Cpj^^^ = 2.16 ± 0.60 mas. The 
O — C residuals from the 5-parameter solution, referred 
to as 'abscissa residuals' Aw, are provided with the cata- 
log and together with the five standard parameters they 
allow the original measurements to be reconstructed. In 
this way, these Aw measurements can be used in principle 
for further improvements in the overall astrometric solu- 
tion if a better visual orbit for Capella were to become 
available. In practice they contribute relatively little to 
the orbit of Capella, but they do provide a useful check on 
the secondary velocity amplitude, to be discussed later. 
Furthermore, they allow an independent estimate of the 
brightness ratio (§131), so we have incorporated these mea- 
surements into our global solution described in the next 
section. They are listed in Table |7l 

Fina lly, Capella was spatially resolved by direct imag- 
ing by lYoung fc Dupreg (|2002|) . using the Faint Object 
Camera aboard HST at ultraviolet wavelengths (1300- 
3000 A) . These measurements are included in Table El 

In many of the interferometric and speckle observations 
the quadrant of the position angles has an ambiguity of 
±180° due to the nature of the measurement. Even in 
cases where the analysis is able to establish the correct 
quadrant, that determination is made more difhcult for 
Capella because the stars are so nearly equal in bright- 
ness, as we discuss in §|5l and because the brightness 
ratio depends on the wavelength of the observation and 
reverses around 7000 A. Here we have adjusted the angles 
where necessary to be consistent with the usual conven- 
tion for visual binaries, in which the position angles are 
measured from the brighter star to the fainter one in the 
V hand. 

Although many of the above astrometric measurements 
have been used previously by others to model the orbit 
of Capella, careful examination during the present work 
of the original references and other bibliographic sources 
making use of them revealed a number of inconsisten- 
cies, misprints, or mistakes that appear not to have been 
noticed before. As a result, the data used here differ 
slightly from a listing of the measureme nts contained in 
the Washington Double Star Catalog (M ason et al.ll200"l[) 
provided by the U.S. Naval Observatory. We document 
these details in Appendix |B1 for the benefit of future 
users. 

4. ORBITAL SOLUTION 

The many data sets described above constrain the pa- 
rameters of Capella's orbit in different ways. While it is 
true that in th i s case the interferometric observations by 
iHummel et al.l (|1994f) and our own radial-velocity mea- 
surements carry much more weight than other data sets, 
the optimal procedure for obtaining the orbital param- 
eters is usually to account for the different weights and 
combine all observations into a single fit, provided they 
are sufficiently free from systematic errors. This is the 
approach we adopt here. The observations consist of po- 
sition angles {9) and angular separations (p), measures 
of the relative separation in rectangular coordinates (Ax 
and Ay) , radial velocities for the primary and secondary, 
and the Hipparcos measurements Av. We solve for the 
usual orbital elements of a visual-spectroscopic binary. 



which are the orbital period (P), relative angular semi- 
major axis (a"), inclination angle (i), eccentricity (e), 
longitude of periastron of the secondary (w), position 
angle of the ascending node for the equinox J2000.0 (il), 
time of periastron passage (T), center-of-mass velocity 
(7), and the velocity semi-amplitudes for each star {Ka 
and Kb)- 

The use of the Hipparcos measurements introduces sev- 
eral additional parameters that must also be solved for. 
These are the angular semimajor axis of the photocenter 
(apjjQj), corrections to the catalog values of the position 
of the barycenter (Aa*, A6) at the mean catalog refer- 
ence epoch of 1991.25, corrections to the proper motion 
components (A/i*, Afis), and a correction to the Hip- 
parcos parallax.^ In this case, however, the fact that the 
spectroscopic elements and Kb are obtained in the 
same solution introduces a redundancy, and the paral- 
lax (referred to here as the "orbital" parallax) can be 
expressed in terms of other elements as 



TTorb = 1-0879 X 10'' 



P{Ka + Kb)VT 



(1) 



The numerical constant is such that the result is in 
the same units as a" (typically mas) when the period 
is given in days and Ka and Kb in kms~^. We have 
therefore chosen to eliminate the parallax correction as 
an adjustable parameter in the fit. The mathematical 
formalism for modeling the Hipp arcos abscissa residuals 
follows closely th a t desc ribed bv Ivan Leeuwen fc Evand 
(1998i), iPourijaixl (|200l . and iJancart et al.l (|2005D . in- 
cluding the correlations between measurements from the 
two independent data reduction consorti a that processed 
the o riginal Hipparcos observations (see iPerrvman et al.l 
Il997f ). Full details along with another exa mple of the ap - 
phcation of this technique may be found in lTorresI (|2007D . 



As noted earlier (i) l2.2p . instrumental effects in spec- 
troscopy often cause the zero points of the radial veloc- 
ity measurements to be different for different observers. 
These shifts are accounted for here by solving for an ad- 
ditional offset between each of the historical RV data 
sets and our own, which we take as the reference be- 
cause it is the largest. We solve for these offsets A, 
(i = 1,...,9) in the sense (other minus CfA) simulta- 
neously with the orbital elements. Additionally, we solve 
for a primary /secondary offset Aab for the CfA veloci- 
ties themselves, to correct for the small shift described 
in ji l2.1l Finally, one more adjustable parameter fp is in- 
cluded as a correction to the scale of the angular separa- 
tion measurements of Merrill (1923) andjKulagin (1970), 
to be described below. Position angles have been pre- 
cessed from the original epoch of e ach observation to the 
standard epoch J2000.0. Those of iHummel et all (| 19941 ) 
have been precessed from their reference epoch of 1991.9. 
For consistency we have also applied precession correc- 
tions to the Aa; and Ay measurements, although they 
are hardly significant. The Hipparcos observations are 
already referred to J2000.0. 

Altogether there are 26 adjustable parameters, which 
we determined simultaneous ly using standard non-linear 
least-squares techniques (see lPress et al.|[l99^ p. 650). A 
total of 1015 individual observations were used in the fit. 

^ Following the practice in the Hipparcos catalog, we define 
Aa* = Aa cos 5 and A/x* = Afict cos 5. 
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A summary of the different data sets can be found in Ta- 
ble [H For approximately half of the astrometric observa- 
tions it was necessary to reverse the quadrants of the po- 
sition angles or the sign of the Ax and Ay measurements 
for consistency. This is hardly surprising given the small 
magnitude difference between the components at optical 
wavelengths, and the inherent ambiguities in quadrant 
determination in some cases. We discuss this further in 
§[5] Uncertainties for the astrometric observations were 
adopted from the original sources, when available, and 
relative weights within each series were accounted for, 
if reported. For some of the speckle measurements that 
have no published errors we adopted typical values of 
fJe — 2° and ap — 3 mas. With few exceptions historical 
radial velocities have no published errors. In those cases 
we have assumed them to be equal to the RMS scat- 
ter from preliminary orbital fits. Relative weights for 
the RVs within a given series were taken into account in 
cases where they were given. Because internal uncertain- 
ties are often underestimated, and some of our guesses 
are necessarily rough, we have re-scaled them by itera- 
tions in the final solution so as to achieve reduced 
values near unity separately for each source, for all as- 
trometric and spectroscopic data sets having a sufficient 
number of observations. 

All prior studies of Capella based on data sets of suf- 
ficient size and quality have concluded that the eccen- 
tricity of the orbit is not significantly different from zero. 
We were therefore somewhat surprised that our initial 
solutions gave a very small yet statistically significant 
value of e = 0.00087 ± 0.00021, with lu = 324° ± 14°. 
Closer examination reve aled that this i s driv en exclu- 
sively by the high- weight iHummel et al.l () 19941 ) observa- 
tions, which when used alone give e = 0.00083 ± 0.00005 
and uj = 334?8 ± 4?7. A solution without the Hum- 
mel measurements yields a circular orbit, as does one 
that uses only the CfA radial velocities, which carry the 
largest weight among the remaining data sets. The CfA 
primary velocities, when considered separately, also sug- 
gest the orbit is circular, but our secondary velocities, 
which have larger uncertainties, prefer e = 0.018 ±0.004. 
This result is clearly related to the residual patterns 
shown in Figure [21 seen only in the secondary, which 
we believe to be most likely of instrumental origin, as 
discussed in ? I2.1I On the basis of this evidence we are 
inclined to conclude that the eccentricity we derive from 
the Hummel et al. measurements is spurious. In their 
own orbital solution those authors made direct use of the 
interferometric visibilities {V^) from the Mark III instru- 
ment, rather than relative positions in polar coordinates, 
which are the data finally published. Nightly values for 
the latter, condensed from the mea sures account- 
ing fo r orbital motion, were provided by iHummel et al.l 
p99l for the convenience of the reader since they are 
easier to use. Given that Hummel et al. reported detect- 
ing no significant eccentricity (e = 0.0000 ± 0.0002) in 
their solutions using the visibilities, we speculate that our 
result is due to our use of the published {d, p} measure- 
ments as opposed to the original values. The transla- 
tion from one to the other has apparently introduced very 
subtle distortions in the orbit, perhaps related to surface 
feature inhomogeneities (spo ts) or calibratio r i issu es, as 
discussed in some detail by IHummel et al.l (|1994[ ). In 
practical terms, the difference between our eccentric and 
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Fig. 5. — Difference between an eccentric and circular or- 
bital solution for Capella, using all data sets. The eccentricity 
is most likely spurio us, and is driven by the measurements of 
IHummel et al.l II1994I ) (see text) . Phases are counted from the time 
of nodal passage in the circular orbit. 



circular fits using all data sets is very small, as illustrated 
in Figure IHl The maximum differences are ~0?1 in po- 
sition angle and ^0.1 mas in angular separation. The 
effect on the absolute masses is considerably less than 
their uncertainties (< 0.5%). For the remainder of this 
paper we will consider the orbit to be circular. This re- 
duces the number of adjustable parameters to 24. The 
epoch T defined above then refers to the nodal passage 
(ascending node) rather than periastron. 

Preliminary fits showed a systematic pattern in the 
residuals of the interferom etric angular separation mea- 
surements of lMerrill (|1922f ). The same pattern is e vident 
in t he orbital solutions published by iMcAfisteH ()1981f ) 
and'Barlo w et al.l (|1993f ) , which show predominantly neg- 
ative residuals in p from this source. Hummel et al] 
(|1994f ) noted a systematic difference between their semi- 
major axis for Capella's orbit and al l previous results , 
beginning with the original study bv lAndersonI l|1920D . 
They speculated that those early interferometric mea- 
surements have a scale problem, and that the large 
weight they have typically received in ot her studies may 
have biased previous orbital solutions. IHummel et al.l 
() 19941 ) also provided a likely explanation for t he scal- 
ing problem. It has to do with the adoption by iMerrilll 
0922) of 5500 A as the effective wavelength used for the 
original Mount Wilson observations. This adopted wave- 
length sets the scale of the angular separations. They 
pointed out that while 5500 A may be a suitable value 
for observations of early G-type stars like the Sun, the 
mean temperature of Capella is now known to somewhat 
cooler than the Sun's, and therefore a slightly longer ef- 
fective wavelength would be more appropriate. In their 
estimation, the early interferometric observations should 
be ^^5% too small. An identical effective wavelength was 
adopt ed in the interferometric observations of Kulagi3 
(|1970f ). and in fact those measurements display the same 
pattern of nega tive r esiduals in the orb ital studies of 
iMcAlisteij (|1981[ ) and lBarlow et all (|1993f ). as weU as in 
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our own preliminary fits. In order to corre ct for this bias 
in t he angular sepa ration measurements of lMerrill (|1922l ) 
and lKulaginI ()1970D . we have included the scale factor fp 
mentioned earlier as an additional free parameter in our 
global solution. Effectively, this means that those mea- 
surements no longer contribute to set the scale of the or- 
bit, but they still help to constrain the remaining orbital 
elements. The result we obtain, fp = 1.0400 ± 0.0035, 
confirms the significance of the effect, whic h is ne arly of 
the magnitude predicted bv lHummel et all |l99^. 

In Table [9] we present our orbital solution for Capella. 
In addition to the adjusted elements, we list a number of 
other properties including the absolute masses and the 
orbital parallax, inferred from the orbital elements. The 
uncertainties for these derived quantities include the con- 
tribution from the off-diagonal terms of the covariance 
matrix, to account for correlations among the elements. 
The determination of the orbital period has benefited 
from the century-long baseline afforded by the observa- 
tions, and its precision is now 2 parts per million (corre- 
sponding to 19.2 seconds out of 104 days). The orbital 
parallax we obtain, TTorb = 76.67±0.17 mas, is consistent 
with, but about 5 times more precise than the value from 
Hipparcos (tthip = 77.29 ± 0.89 mas).^ 

Residuals from the spectroscopic observations are pre- 
sented in Table [T] and Table [3l while those of the astro- 
metric observations are given in TablelH TableO Tabled 
and Table [71 The typical precision of the measurements 
from each source as represented by the RMS residual of 
unit weight is given in Table [H As indicated earlier, the 
Mark III observations by Hummel ct al. (1994) are by far 
the most precise of the astrometric data, and are shown 
graphically in Figure [6] separately from the other obser- 
vations. Residuals in position angle and separation are 
also shown, and are typically 0?13 in 9 and 0.11 mas in p. 
The speckle observations are displayed in Figure [71 with 
their residuals shown on the same scale as the previous 
figure for comparison. All other measurements obtained 
by long-baseline interferometry are given in Figure [H in- 
cluding both those originally made in polar coordinates 
and those made in rectangular coordinates. Among the 
latter, the much larger residuals in Aa: (right ascension) 
than in Ay (declination) are due to the north-south ori- 
entation of the basel ine o f the interferometer used by 
iKoechlin et"all (|1979D and lKoechhn et all (|1983D . which 
is the source of most of those measurements. 

Examination of Table [H r evea l s that the residuals in 9 
for the observations by iMerrilll ()1922l ) show a tendency 
toward negative value s for the later d ates. The earlier 
ob servations ma de bv [Anderson! (|1920( ) (and re-reduced 
bv [Merrill 1 1 92 2f l show the opposite trend, with the ex- 
ception of the very first measurement, which is of much 
lower quality and has a very large error. These trends 
were noticed already by Merrill, who offered as explana- 
tions either an instrumental effect or a real advance of 
the node. We find no evidence for a secular change in 
J7 in the other observations, so we tend to agree with 

^ A recen t new reduction of the Hipparcos observations by 
Ivan LeeuwenI I I2007I) yielded an improved parallax value for Capella 
of TTHip = 76.19 ± 0.47 mas, which was subsequently revised in 
the online version of the catalog to correct for an error that af- 
fected the goodness of fit in some cases. The updated value, 
TTHip = 76.20 ± 0.46 mas, is still within Itr of our more precise 
determination. 




1 . , . . 1 . , . . 1 . 


50 


50 
Aa (mas) 


- 1 ' 1 1 

- 1 , 1 


1 ' 1 ' 1 1 1 - 
• •• --•--m «» • 

1 , 1 , 1 , 1 - 




_ 1 ' 1 

" 1 , 1 


1 ' 1 ' 1 ' 1 _ 
• •• • «• mm «• A • 

1 , 1 , 1 , 1 


0.0 0.2 


0.4 0.6 O.B 1.0 
Orbital phase 



Fig. 6. — Mark II I interferometric observations of Capella by 
Humm el et al.l I I1994I) . together with the orbit computed from our 
global solution. The visually brighter component (star B) is at 
the center. Solid lines connect the observations with the predicted 
position on the orbit. O — C residuals in position angle and sepa- 
ration are indicated at the bottom. Error bars are mostly smaller 
than the size of the points. The line of nodes is indicated with the 
dashed line, and motion on the sky is retrograde (clockwise). 



[McAlisteil ([1981[ ) that it is most likely an instrumental 
problem.'' As a test, we repeated the orbital solution 
solving for two position angle corrections in addition to 
the other 24 elements. We obtained — 0?5 ± 0?5 for the 
earlier observations by Anderson, and -t-l?4±0?4 for the 
later ones by Merrill, consistent with expectations. Ad- 
justing the original values of 9 for these offsets leads to 
a la decrease in the orbital period of Capella, and a 
slightly reduced uncertainty in P of 17 sec. The change 
in all other elements and derived quantities is negligible. 

Given that the components of Capella are slightly dif- 
ferent, the size of the apparent orbit described by the 
center of light of the binary as seen in unresolved obser- 
vations depends on the wavelength of the observation. 
Our inclusion of the Hipparcos data in the solution en- 
ables us to derive the brightness ratio Ib I^a between the 
stars in the passband of the satellite, denoted Hp. For 
this we make use of the fact that the semimajor axis of 



^ IMerrilll II1922I ) himself pointed out that there was no direct 
way of checking the position angle circle of the instrument when 
attached to the telescope, so that the actual position angles of the 
interferometer slits could have differed by small amounts from the 
angles as read from the circle. 
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Fig. 7. — Same as Fig. [6] for all speckle observations of Capella. 
For reference, the orbit of the center of light of of the binary as 
seen by Hipparcos is shown to scale in the upper right corner (see 
also Fig.[9]|. O — C residuals are indicated at the bottom on the 
same scale as in Fig.|6l to facilitate the comparison. 

the photocenter and that of the relative orbit are related 
by aj^'hot = - where B = Mb /(Ma + Mb) is the 
fractional mass and B = £b/(^a + ^b) i s the fractional 
luminosity (see, e.g.. Ivan de Kamdfl967l ). This leads to 



{iB/e_ 



A)H„ - 



Ka 



^phot 



Ka + Kb 



- 1 



(2) 



Our resulting light ratio along with other estimates of the 
relative brightness are discussed in §[5] The projection 
of the photocentric orbit of Capella on the plane of the 
sky along with a schematic representation the Hipparcos 
measurements is seen in Figure [SI The much smaller size 
of the photocentric orbit compared to the relative orbit 
is illustrated in Figure [71 

5. THE LIGHT RATIO 

The near equal brightness of the components of Capella 
has been a source of considerable confusion in the 
past. The quadrant of the ascending node and the time 
of nodal passage (or equivalently, the identity of the 
brighter star) have been changed more than once since 
the public ation of the first set of astrometric orbital el- 
ements by lAndersonI ()1920D .^ The spectrophotomctric 

* The choice of q uadrant in th at work appears, however, to have 
been arbitrary (see IFinsenlll975l ) . 
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Fig. 8. — Same as Fig. [6] and Fig. [7] for all long-basclin e inter- 
ferom etric observations of Capella except those of Humm el et al.l 
l|T99l ). O-C residuals are shown at the bottom, separately for 
observations made originally in polar coordinates {9, p ; filled cir- 
cles) and in rectangular coordinates (Ax, Ay ; open circles). The 
scale of the position angle and separation residuals is the same as 
in Fig. [6] and Fig. [7] to facilitate the comparison. 



study bv IWrighB (|1954f ). in which the author incorrectly 
concluded that the cooler star was the brighter one in 
the visible by ^0.25 mag, played an important role in 
our understanding of the system for several decades, al- 
though unfortunately it also introduced biases in a num- 
ber of other investigations that made use of that re- 
sult. Exampl es include, among ot hers, the interferomet- 
ric study by iBlazit et al.l (|1977a[ ). who attempted the 
first angular diameter measuremen ts of the stars, the 
Li abundance de terminations by iWallersteinI fl966. ) and 
iBoesgaardI (|1971[ ). and to some extent the ^^C/^^C ratio 
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Fig. 9. — Motion of the photocenter of Capella relative to the 
center of mass of the binary (indicated by the plus sign) as seen by 
Hipparcos . See Fig.[7]for a comparison with the size of the relative 
orbit. The solid curve is the computed orbit from our global solu- 
tion. Because these measurements are one-dimensional in nature, 
their exact location on the plane of the sky cannot be displayed 
graphically. The abscissa residuals are shown schematically with a 
filled circle at the predicted location, dotted lines representing the 
scanning direction of the satellite, and short perpendicular line seg- 
ments indicating the undetermined location of the measurement on 
that line. The length of the dotted lines represents the magnitude 
of the O — C residual from the computed location. Measurements 
with large residuals have been omitted for clarity. Motion is retro- 
grade (clockwise). 



estimate by iTomkin at al.l ()1976f) . all of which adopted 
Wright's brightness ratio. The history of this problem 
has been well summarize d by | Grif5n fc Griffi.nl (|1986D . 
and further discussed by iBarlow et al.1 ( 19931 ) . Begin- 
ning in the early 1980s a number of authors used long- 
baseline interferometry and speckle interferometry tech- 
niques to unambiguously identify the hotter star as the 
brighter one in t he yis ible (shortward of 7000 A), and 
I Griffin fc GriffinI (|1986f ) provided a reasonable explana- 
tion for Wright's spectroscopic result, which apparently 
referred to a difference between continuum heights rather 
than relatiye light intensities, and did not account for the 
difference in line blocking between the stars. 

Nearly two dozen independent measurements of the 
relative brightness of the components are now available 
from the ultraviolet to the infrared. They are listed in 
Table [TOl and are combined later with absolute pho- 
tometry to derive effective temperatures for the compo- 
nents. Included also are our own light-ratio estimates 
from the CfA spectroscopy (see § 12. ip and from the Hip- 
parcos observation s (§B1). The spectroscopic measure- 
ment reported by IStrassmeier fc Fekel (|1990[ ) refers to 
the ratio of continuum heights, rather than the inten- 
sity ratio. To convert to a true intensity ratio we have 
applied a correction for the line blocking based on ap- 
propriate synthetic spectra over the passband of their 
observations. We report the corrected value in Table [TOl 
For uniformity the quantities listed in the table have all 
been converted to ratios ^cooi/^hot {— ^a/^b) between the 
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Fig. 10. — Light ratio between the components of Capella as a 
function of wavelength. The shaded regions represent the uncer- 
tainty (vertical direction) and the FWHM of the passband (hori- 
zontal direction). 

cooler, slightly more massive star and the hotter star. A 
graphical representation of these measurements is seen 
in Figure [ini At wavelengths near the V band the hot- 
ter star is slightly brighter, but around 7000 A the ra- 
tio reverses, and the cooler, more massive star becomes 
dominant. This explains why some of the interferomet - 
ric measurement s such as those by IB aid win et alj ()1996f ) 
and lKraus et al.l (|2005D . which were made at red or near 
infrared wavelengths, have the quadrants reversed com- 
pared with measurements in the optical. 

6. ANGULAR DIAMETERS 

The angular sizes of the components of Capella are 
large enough that they have been resolved by long- 
baseline interferom etry on several occas ions. They were 
first measured by iBlazit et al.l (|1977al) . who obtained 
uniform-disk angular diameters of 9^^ = 5.2 ± 1.0 mas 
for the cooler primary star and — 4.0 ± 2.0 mas for 
the hotter star. However, these values assumed that the 
co oler star is bri ghter by 0.25 mag, following the results 
of I Wrightl (|1954l ). whereas we now know the cooler star 
is in fact the fainter one (see §[5]). Unfortunately it is 
not possible to correct Blazit' s original estimates based 
on the information reported. iDi Benedetto fc Bonneaiil 
(|1991f ) obtained a limb-darkened angular diameter of 
®LD ~ 6-2^ ^ 0.43 mas for the secondary, along with 
a much more uncertain value of Old ^ 9.6 ± 2.3 mas 
for the primary, both measured in the H and K bands. 
Un iform-disk di a meter s in the H band were published 
by iKraus et al.1 (|2005f) as = 8.9 d z 0.6 mas and 

OPin = 5.8 ±0.8 mas. The observations of lHummel et al.l 
(|1994l ) at wavelengths corresponding approximately to 
the -B, V, and Iq bands gave limb-darkened angular 
diameters for Capella of Old = 8.5 ± 0.1 mas and 
O^L) = 6.4 ± 0.3 mas. 

The above measurements are inhomogeneous due to 
the variety of limb -darkening correc t ions u sed. Those 
apphed by ,Di Benedetto fc Bonneaul (|199lD correspond 
to a scale factor of 1.035 between Old and Oud. 
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following Carpented ()2001f ). we obtain 
mas and 



iHummel et al.l (119941) used limb -darkening coeffi cients 
from iManduca et all ( 1 19771 ) and iManducaT ()1979f ). and 
iKraus et all (|2005D chose not to apply any corrections 
at all. To place all these measures on the same foot- 
ing we have a dopted limb-d arkening coefficients from the 
tabulation by Ivan Hammd lil993.1. and cornputed the 8 
corrections following iHanburv Brown et all ()1974[ ). The 
differences in these corrections compared to the origi- 
nal ones can be as large as 1.7%. The homogenized 
angular diameters are listed in Table [TlJ The result- 
ing weighted averages are = 8.47 ± 0.40 mas and 
^LD " 6-24 ± 0.23 mas. The uncertainties, which ac- 
count for the scatter in the individual measurements, 
correspond to fractional errors of 4.7% for the primary 
and 3.7% for the secondary. These angular diameters, 
combined with the orbital parallax, yield the absolute 
radii of the components that are presented below. 

As a check, independent estimates of the angular di- 
ameters may be obta ined from the r iear-i nfrared surface- 
brightness relation of Di Be nedettol (|1998D for giant stars, 
which is very tight and has a scatter of only 1.4%. The 
required V—K indices for the components of Capella are 
available from published photometry and are described 
in the next section (see also Table fT3|). After transforma- 
tion of th e photome t ry to the standard Johnson system 

8.32±0.75 

6.17 ±0.65 mas, in which the uncertain- 
ties include all photometric errors as well as the scatter 
of the calibration. While less precise, these values are in 
excellent agreement with the direct measurements from 
interferometry. 

7. CHEMICAL COMPOSITION 

Chemical composition plays a very important role in 
the comparison with models in the following sections, 
and provides important clues on the evolutionary state 
of the system. In this section we critically review and 
discuss all available abundance determinations in some 
detail, most of which have never been used before in the 
analysis of this binary. 

Despite being such a bright star, the determination 
of the photospheric chemical composition of Capella 
has received relatively little attention by spectroscopists. 
The on ly detailed high-re solution study appears to be 
that of iMcWilligiml (|1990f) . in the context of a survey 
of 671 G and K giants. The value reported is [Fc/H] 
= -0.37 ±0.22 on the scale of lGrevessd (|1984D . in which 
the abundance of iron is logiV(Fe) = 7.67. This re- 
sult is presumably based on the sharp lines of the pri- 
mary. It does not seem that the study has accounted 
for the double-lined nature of the spectrum, which can 
influence the metallicity significantly in two ways. On 
the one hand, the continuum of the secondary (which 
has the same brightness a s the o ther star at the wave- 
lengths of the iMcWiUiaml ()1990f ) spectra) wiU tend to 
fill in the lines of the primary at most phases, making 
them look weaker. On the other hand, the temperature 
adopted for Capella in this analysis (5270 K) was based 
on the combined-light photometry, and is too hot if as- 
signed solely to the primary. This will generally result 
in abundances that are too high. It may compensate 
for the other effect to a certain degree, but the net bias 
is difficult to predict. Aside from the particular case of 
Capella, small systematic differences in the iron abun- 



dances between this w ork and others have occasionally 
been pointed out (e.g.,lLuck fc WepfeHll995t I Zhao et al] 
120011: Ida Silva et al.ll2006l: iLiu et al.ll2007[ ). and are prob- 
ably traceable to systematic differences in the adopted 
surface gravities or microturbulent velo cities. 

Here we place the iMcWilliamI ()199Q) [Fe/H] determi- 
nation on the more recent scale of solar abundances by 
iGrevesse fc Sauval (|1998D . used in some of the models 
considered later, in which logiV(Fe) = 7.50. We ob- 
tain [Fc/H ] = —0 . 20 ± 0.22, where the error is repeated 
from McWi lliamI (jl990r ) and corresponds to the scatter 
of the individual iron line measurements rather than the 
uncertainty of the mean. The abunda.nces of a dozen 
other elements studied by iMcWilliamI ()1990r ) were sim- 
ilarly converted to the same scale, and are collected in 
Table 1121 For the purpose of comparison with the mod- 
els in j; l9.1[ which assume sol ar-scaled abundances, we 
follow I Valenti fc Fischeil ()2005D and adopt the average of 
all elements as an overall indicator of metallicity: [m/H] 
= —0.34 ± 0.07. The uncertainty given here is the error 
of the mean. There is no evidence for enhancement of 
the a elements. All other indicators of the photospheric 
composition of Capella found in the literature are either 
circumstantial, contradictory, or inconclusive.^ 

The photospheric ^^C/^'^C isotope ratio has been mea- 
sured in the optical for the primary star bv lTomkin et al] 
(|1976f ). who reported the value 27 ± 4. A s noted e arlier 
this study used the incorrect light ratio of lWrightl ()1954D 
to subtract the contribution of the secondary to the ob- 
served continuum, although we do not believe this intro- 
duces a large error due to the differential nature of the 
measurement. "'^^ This isotope ratio is a valuable indicator 



^ lEggenl | |T960| . \l97^ regarded Capella as a member of the 
Hyades moving group, primarily based on kinematic criteria. We 
confirm that assessment: we obtain UVW velocities of —36.5, 
— 13.9, and —9.1 kms"'^ (with U toward the Galactic center), in 
good agreement with the mean values and dispersions f or the group 
of -3 8 ± 6, -17 ± 6, and -11 ± 12 kms'^ given bv IZhao et atl 
II2009I) . This circumstantial evidence would imply a composition 
near solar, since the mean metallicity of the group ap pears to be 
[Fe/H] = -0.09 with a scatter of 0.17 dex i fZhao et atll2009t) . Un- 
fortunately our own spectroscopic material does not allow an accu- 
rate determination of [m/H] because of the strong dependence of 
metallicity on temperature over the narrow wavelength range avail- 
able (see §[2]l- Other estimates of the photospheric abundance scat- 
tered throughout the literature show very poor agreement. A rough 
determination by Miner ( 1963) based on photometry using narrow- 
band interference filters ga ve an overa ll composition near solar for 
the combined light. Boesga ardI 1 119711) measured the Li abundance 
of Capella, and in the same study listed also an iron abundance 
of [Fc/H] = 4-0.26. Few details of this determination were given, 
aside from the fact that the equivalent widths of the iron lines for 
each component were corrected for the ligh t contribution from the 
other star using the light ratio of IWrigM 0-954), which we now 
kn ow to be r eversed (see §[5]l. In their lithium study of Capella 
Pilachowski & Sowoll ( 1992) did not report an iron abundance, but 
pointed out that t he calcium abundance is essentially solar for both 
components. Ra ndich et al.l II1994D reported [Fe/H] = —0.4 for the 
primary of Capella, and solar metallicity for the secondary. They 
speculated the discrepancy could be due to differences in chromo- 
spheric activity, although they also noted that other evidence goes 
against this. Fi nally, a s tudy of the coronal metallicities from X-ray 
observations bv lBaueT^ Bregman (1996) mentions a photospheric 
metallicity corresponding to [Fc/H] = +0.27, and attributes this 
determination to Mercki, Strobel & Strobel (1986) without giving 
a bibliographic reference. We are unable to trace this source in the 
literature. 

We estimate the equi valent width measurements of the CN 
lines reported by Tomkin e t al.l II1976I) to be underestimated by 
about 6—11% due to this efffect. 



Capella 



13 



of evolution. 

A large difference in the lithium abundance be- 
tween Cap e lla A and B was first pointed out by 
IWallersteirJ ()1964[ ). and confirmed by others. The hot- 
ter secondary has approximately two orders of mag- 
nitude stronger lithium than the primary. Measure- 
ments h ave been made by|W allcrstcin (1966), Boesgaard 
(119711). iPilachowski fc Sowell (,1992, ). ,Liu et al., (,19931 ). 
and iRandich et aL I (Il994f) . in which the first two are af- 
fected by the use of Wright's light ratio, and the latter 
three adopted effective temperatur es somewhat different 
from o urs. The measurements by IPilachowski fc Sowelll 
()1992D appear to be the most reliable, although the oth- 
ers are generally consistent when a djusted for the modern 
light ratio. Here we have used the IPilachowski fc Sowelll 
(fl992. ) equivalent width measurements for the Li I 6708 A 
feature (25 ± 2 mA and 200 ± 10 mA for the pri- 
mary and secondary, respectively) . We recalculated the 
abundances using the models by iPavlenko fc Magazzul 
()19 96^. accounting for the temperature and gravity dif- 
ferences as well as non-LTE effects (not considered in 
the original analysis). We obtain revised lithium abun- 
dances of logA^(Li) = 1.3 ± 0.2 for the primary and 
logA^(Li) = 3.2 ± 0.3 for the secondary, in which the 
uncertainties include all measurement errors as well as 
possible errors in the microturbulent velocity following 
IPilachowski fc Sowelll ((l992). To be conservative, the 
uncertainties have been further increased by 0.1 dex to 
account f or slight extrapolations that were necessary in 
using the IPavlenko fc Magazzijj |1996) tables. 

As an active binary system, Capella has been stud- 
ied extensively in the ultraviolet and X rays for decades 
using virtually every space facility capable of observ- 
ing at those wavelengths. It was in fact the first X- 
ray detection of a stellar coro na other than the Sun, 
made by sound ing rockets ( Catura et al.l [l975l : see also 
iFisher fc Mever ott 1964; Avrcs et al. 19 91). At u ltravi - 
olet wavelengths, .Bohm-Vitense fc Mena-WerthI (|1992D 
have presented evidence that reliable abundance ratios 
between carbon and nitrogen can be determined for gi- 
ant stars from measurements of the emission fluxes of 
the C IV A1550.8 and N V A1238.8 lines in the lower 
transition layers between the stellar chromosphere and 
the corona, and that these ratios show good correspon- 
dence with the photospheric abundance ratios. Emission 
fluxes for these lines have been measured in Capella by a 
number of authors. However, early observations did not 
clearly resolve the contribution of the two components, of 
which the p rimary represents on ly ^^10%. This was first 
achieved by iLinskv et al.l (jl995') based on high spectral 
resolution observations with HST. Using the fluxes they 
reported, we have derived the C/N ratios for Capella and 
use them below as diagnostics of evolution. 

Abundance determinations for Capella have also been 
made by many authors from X-ray observations of coro- 
nal lines. With current instrumentation it is gener- 
ally not possible to separate the spectral contribution 
of the two componen ts in X rays, as it is in the UV. 
Ilshibashi et al.l (|2006D and others have reported that the 
cool primary often dominates the coronal emission in 
this spectral region, although its flux is variable with 
time. Others find a more nearly equal contribution (e.g., 
ILinskv et aLlll998f) . Therefore, any measurements may 



refer mainly to the primary, but are likely to be contami- 
nated by the secondary. Most of these observations have 
revealed e nhanced nit rogen (see, e.g., 'Mewe et al.' '2001|: 
I Schmitt fc Ness 200l lAudard et al.l l2003: Argiroffi et aT| 
|2003[ ). Mean abundances averaged over all other ele- 
ments from these studies typically indicate a sub-solar 
composition, in qualitative agreement with t he pho- 
tospheric determinations of iMc William! ()1990[ ). How- 
ever, coronal metallicit y measurements are hi ghly model- 
dependent (see, e.g.. iBrickhouse et al.|[200ll ). and indi- 
vidual values sometimes show a large scatter from au- 
thor to author. Furthermore, in the Sun's coronal re- 
gions abundances are known to depend on the first ion- 
ization potential (FIP) o f the element considered (e.g., 
iFeldman fc Widind l200l . In view of these complica- 
tions, we have preferred not to make use of these data 
here. Nevertheless, the nitrogen enhancement consti- 
tutes an interesting piece of chemical evidence for the 
evolved state of the primary, as recognized by many au- 
thors, since it is a natural consequence of the CNO cy- 
cle for stars that have already experienced first dredge- 
up (see below). Careful consideration of the FIP effect 
in Capella suggests there may be even closer agreement 
between the overall coronal abundance and the photo- 
spheric value, which we believe is worth noting given our 
concerns expressed earlier about the latter. We discuss 
these coronal measurements and their patterns in Ap- 
pendix [C] All other useful abundance determinations 
for Capella described above are gathered in Table [121 

8. ABSOLUTE DIMENSIONS 

The astrometric-spectroscopic orbital solution in §[4] 
yields directly the absolute masses of the components. 
The relative uncertainties (0.7% and 0.5 %) represent a 
factor of 3 improvement over those of iHummel et al.l 
(|1994| ). which is critical for the comparison with stellar 
evolution models. We also obtained the orbital paral- 
lax. The formal uncertainty in the corresponding dis- 
tance of 13.042 ± 0.028 pc is only 0.2%. The absolute 
radii of the components follow from the angular diame- 
ters and the distance, and are Ra = 11.87±0.56 Rq and 
i?B = 8.75 ±0.32 Rq. Relative to the orbital separation, 
these values correspond to 0.075 and 0.055, respectively, 
so the binary is well detached. 

Effective temperatures for the individual stars in 
Capella have been estimated here in three different ways. 
A first determination relies on our spectroscopy, and is 



described in §0 



4900 ± 100 K and T^l 



5710 lb 100 K. A sec ond method is that employed by 
iHummel et al.l ()1994[ ). who made use of their angular 
diameter measurements along with the apparent mag- 
nitudes and bolometric corrections to infer values of 
r_A ^ 4940 ± 50 K and Tjf = 5700 ± 100 K, nearly 
identical to ours. We have updated that calculation 
using the average angular diameters from §|6l together 
with apparent visual magnitudes for the components as 
described below, and bolometric corrections BCv from 
Flower (1996)." The results are = 4970± 154 K and 
= 5687 ± 130 K, in which the uncertainties in BCy 

To be consistent with the scale of the bolometric corrections, 
the bolometric magnitude adopted here for the Sun is Af^j = 
4.732. When combined with the tabulated BCy corresponding to 
the solar temperature of T^g = 5777 K, this gives an apparent 
magnitude for the Sun that reproduces the measured value of y = 
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and in all other measured quantities are included. A 
third method to estimate individual temperatures relies 
exclusively on photometry (color indices), and has been 
applied by a number of authors over the years giving re- 
sults generally consistent with the above estimates. We 
return to this technique below. To our knowledge there 
are no other fundamentally different Tes estimates avail- 
able, except for those one might infer indirectly from the 
spectral types assigned to t he components. For example, 
IStrassmeier fc Fekell ()199 0') applied a spectrum synthesis 
technique and found a good match to the primary and 
secondary in the standard stars Pollux {(3 Gem, KO III) 
and a Sge (Gl III). These classifications are roughly con- 
sistent with our Teg estimates. 

The information on the absolute photometry for 
Capella and the light ratios discussed earlier is col- 
lected in Table [T31 We use it here to derive photo- 
metric estimates of Tcs for each star. The light ratios 
in the table are weighted averages of all values near 
the B, V, R, I, J, H, and K passbands, respectively, 
and the co mbined-light magn itudes are taken from the 
database of iMermilliod et al] ([l997). R and / magni- 
tudes w" e re tra nsformed to the Cousins system following 
iLeggettI (|1992D . and the near-infrared magnitudes were 
place d on the 2MASS system using the transformations 
of Ca rpenter! ([2001). Individual uncertainties are taken 
as published. Color indices formed from these values are 
listed in the second section of the table. These are not 
strictly independent, but they at least provide a sense of 
the consistency of the measurements in different systems 
and the scatter one can expect from the external calibra- 
tions applied in each case. Color /temperature ca libra- 
tions for giant stars bv lRamirez fc MelendezI ((2005h were 
used to derive temperatures for each component as well 
as for the combined light (third section of Table [T5)). The 
metallicity adopted is the value [m/H] = —0.34 zfc 0.07 
based on the measurements by IMcWilliamI (|199Clf ). de- 
scribed in the previous section. The temperature uncer- 
tainties reported in the table account for all photometric 
errors, the uncertainty in the assumed [m/H], and also 
the scatter of each color/temperature relation. Weighted 
average temperatures computed from the seven indices 
are listed as well. The values for both Capella A and B 
are in good agreement with the other two determinations 
described previously. 

The last line of Table [13] presents the weighted 
average of the three independent Tcs determinations 
for each star, based on the spectroscopy, the quan- 
tities {Old, V, BCv}, and photometry, respectively. 
To be conservative, the uncertainty of the photomet- 
ric values have been increased by adding 100 K in 
quadrature to the formal errors prior to taking the 
average, in order to account for possible systemat- 
ics in the color / temp erature calibrations. This fol- 
lows the discussions ofjR amirez fc Melended ^mh and 
ICasagrande et al.l (|2006D concerning our knowledge of 
the absolute effective temperature scale. The final tem- 
peratures are 4920±70 K and 5680±70 K for the primary 
and secondary, respectively. 

The very different rotational velocities of the compo- 
nents was already evident to spectroscopic observers a 

-26/76 ±0.02 as determined bv lStebbins fc KronI I^T957^ and lHaved 
II1985I '). See also the discussion bv lBessell et al.1 II1998I') . 



century ago. The vsini values have since been measured 
by many investigators, mostly by traditional spectro- 
scopic means but also with other methods s uch as the dif- 
ferent ial speckle interferometry technique of lPetrov et al.l 
11996') . These estimates are collected in Table [M] along 
with our own. For the most part the more recent deter- 
minations agree fairly well, considering the difficulty of 
the measurements. 

The physical parameters for both components of 
Capella are summarized in Table 1151 The luminosi- 
ties were derived here from the well determine d abso- 
lute m agnitudes and bolometric corrections from lFloweil 
(Il996l). The uncertainties in BCy were propagated from 
the error in Tgff, and an additional conservative error 
of 0.05 mag was added in quadrature. If we instead 
compute the luminosities directly from the radii and 
temperatures, the values are considerably more uncer- 
tain (La = 74.2 ± 8.2 Lq, Lb = 71.5 ± 6.1 Lq), but 
are consistent with the adopted estimates. The primary 
(cooler) star is the more luminous bolometrically, but is 
the fainter one in the visible. Also included in the table 
are the projected rotational velocities (wgyncSini) com- 
puted under the assumption that the stars have their 
spins synchronized with the orbital motion and that the 
spin axes are perpendicular to the orbital plane. We dis- 
cuss these values in §|9l 

9. DISCUSSION 

The key properties that determine the evolutionary 
state of the giants in Capella are the masses. Prior 
to this study the val ues most often adopted (e.g., 
iNobuvuki fc Saiol fl999l) were those of iHummel et al.l 
([199J), Ma = 2.69 ± 0.06 Mq and Mb = 2.56 ± 
0. 04 Mp,, which rely on the velocity semi-amplitudes 
of iBarlow etHl (|1993[ ). These masses are 9% and 5% 
larger, respectively, than those in the present work. As 
noted in § 12.21 our primary velocity semi-amplitude is 
not very different from other determinations, but our 
secondary value is considerably smaller, and this drives 
both masses down. An independent check on the ac- 
curacy of Kb can be made with the available astrom- 
etry (specifically, the Hipparcos observations), without 
using any secondary velocities. This is because the Hip- 
parcos measurements are on an absolute frame of ref- 
erence (ICRS) and therefore contain strong information 
on the trigonometric parallax, and the parallax is related 
to Kb via eq.([T]). We carried out an orbital solution in 
which our secondary velocities were given zero weight, 
and the result for the secondary semi-amplitude is Kb = 
26.01 ±0.62 kms~^. This is considerably more uncertain 
than the spectroscopic value of 26.260 ± 0.087 kms~^, 
but is still perfectly consistent with it, while at the same 
ti me being more than 2a away from the determination 
by iBarlow et al.l ([19931 ) . This suggests our masses for 
Capella are more accurate than previously determined, 
in addition to having smaller formal errors, and we pro- 
ceed below to compare them along with other observa- 
tions against theory. 

9.1. Comparison with stellar evolution models 

Detached binary systems such as Capella that are com- 
posed of two giant stars and show double-lined spec- 
tra are rare, and they provide important tests of mod- 
els in a relatively short-lived phase of stellar evolution. 
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Their component masses are necessarily very close to 
each other, and a precise measurement of the mass ratio 
q, as we provide here, becomes critical to establishing 
their state of evolution unambiguously. 

The evolutionary status of Capella has been a subject 
of debate for decades. While there is general consensus 
that the hotter secondary is crossing the Hertzprung gap 
and approaching the base of the giant branch (RGB), 
opinions have varied on the precise location of the pri- 
mary, in large part because of uncertainties in the masses 
as well as the effective temperatures and luminosities 
used to place the star on the H-R diagram. Capella is 
perhaps unique in that, in addition to those properties, 
a wealth of other information is available to aid in de- 
termining its evolutionary state, including the surface 
Li abundances of both stars, the ^^C/^^C isotope ratio 
for the primary, C/N ratios, and activity indicators in 
the optical, ultraviolet, and X rays. The progenitors of 
Capella were late B- or early A-type stars. When such 
stars leave the main sequence, they develop convective 
envelopes that deepen significantly as they approach the 
giant branch, mixing the outer layers with matter from 
the interior partially processed through the CNO cycle. 
As a result, fragile elements such as lithium are burned 
deeper in the star decreasing the surface abundance of 
that element, and others such as ^'^C and -'^"'N that are 
created at the expense of ^^C are brought to the surface 
during the "first dredge-up" . This causes a dramatic re- 
duction in the ^^C/^^C ratio and in the C/N ratio, both 
of which are measurable. Thus, surface abundances con- 
tain potentially important information on the evolution- 
a ry sta t e of ev olved stars like Capella. 

Ilbeni ()1965l ) pioneered this approach by relying on 
early est imates of the lith i um a bundance of both com- 
ponents (|Wallersteinlll96 j . Il966f) to conclude, based on 
his models, that the primary is a core helium-burning 
star. Compared to an alternate location on the ascend- 
ing giant branch, the "clump" phase also seems more 
likely because it is longer-lived.^^ This lifetime argu- 
ment seems to have weighed heavily in most of the 
other studies in which the measured masses, tempera- 
tures and luminosities have been compare d against stel- 
lar evo l ution models, in c luding t he work of iBarlow et a.Ll 
(I1993D. iHummel et all (Il994li . ISchroder et"al] (|1997f ). 
and [iw amoto &: Said ()1999D . On the other hand, 
iBoesga ard (1971^ concluded based on her own Li mea- 
surements, which differed from those of Wallerstein, that 
the primary is not in such an advanced e volutionary 
state. Similarly, iBagnuolo fc Hartkopj ()1989D found evi- 
dence in the small luminosity difference between the stars 
that the primary is still at the beginning of the RGB. 
They also argued for a much smaller difference in mass 
than indicated by the measurements at the time, and 
indeed our present determinations bring the mass ratio 
much closer to unity than implied by the spectroscopy 

The predicted durations of the diflferont stages of evolution 
according to one of the models of Claret ( 2004) considered below 
(case A) are as follows, for a star with the mass of the primary. 
The main-sequence (MS) phase lasts 526 Myr. The crossing of the 
Hertzprung gap up to the point of minimum luminosity at the base 
of the RGB is 7.3 Myr, or only 1.4% of the MS Ufctime. The first 
ascent up to the helium flash lasts 5.7 Myr (1.1%), the subsequent 
descent to the luminosity minimum takes 16.2 Myr (3.1%), and the 
clump phase is a more prolonged 89.8 Myr (17.1%). 



of IBarlow et"all (|1993D . lAvres et al.l (|1983[ ) also took the 
view that the primary is not yet burning helium in its 
core based on the high levels of chromospheric activity 
implied by their ultraviolet observ ations, although the 
opposite conclusion was r eached bv lAvresI (|1988l ). 

The significantly improved parameters we have derived 
for Capella, particularly the masses which are three times 
more precise than those previously available, offer an 
opportunity to revisit the issue of its evolutionary sta- 
tus. Among the many publicly available stellar evolution 
models, we initially focused on those that allow some flex- 
ibility in setting parameters such as the composition or 
the age (for isochrones). However, not all of these models 
extend past the helium flash, which is necessary to ex- 
plore the more advanced core helium-burning phase for 
the primary, and thus we are somewhat limited in our 
choices. 

As a starting point, we compare the observations 
against the widely used set of stellar models by 
iGirardi et al.l (|2000f ) , based on physics that are now stan- 
dard in most current models, including convective core 
overshooting. The mixing length parameter is fixed in 
these models at the value auh = 1.68, overshooting is 
set to Q!ov = 0.25, and the chemical composition adopted 
for the Sun is Zq — 0.019. Diffusion is not considered, al- 
though its effect is completely negligible for Capella. Nei- 
ther is mass loss due to winds, whi ch appears to be rather 
low according to most estiinates fPrake fc Linskvl [T986I : 
iKatsova fc Scherbakovlll998l : iGetman fc LivshitsI 11999( 1. 
For convenience we have chosen to compare the observa- 
tions against isochrones (thereby imposing the constraint 
of coevality) computed using the web interface provided 
by the authors. In addition to M, L, and Tes, we con- 
sider also the absolute radii since they are determined 
independently from the temperatures and luminosities 
and are of comparable precision. We explored a fine grid 
of ages and chemical compositions, and used as a figure 
of merit the defined as 




in which the A quantities represent the differences be- 
tween the measurements and the models for each star 
(i = 1,2). The best match to the observations is achieved 
for a metallicity Z = 0.008 (corresponding to [Fe/H] = 
—0.38) and an age of r = 537 Myr. All four measured 
quantities are reproduced to within 1.4 times their nomi- 
nal uncertainties (see Table [TBI case A), and the primary 
star is located on the ascending branch. A similarly 
good fit to an isochrone of the same age and composi- 
tion can be found with the primary on the descending 
branch (Table [HI case B), with only a slightly larger 
difference with that star's measured effective tempera- 
ture (l.Tcr).^"' These two fits are shown in Figure [TTl 
No satisfactory match is found with the primary in the 
core helium-burning phase, which has been favored by 

|http: // stev. o apd. Inaf . it/ cgi-bin/cmd| . 

We note that the mass ratio for Capella is better determined 
from our measurements (0.36% relative error) than the individual 
masses (0.73% and 0.53%, respectively). The best fit with the 
primary on the ascending branch predicts a value of q which is 
2.2cr larger than measured, while that on the descending branch 
shows better agreement at the 0.3iT level. 
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Fig. 11. — Isochrone from the series bv lGirardi eFall l[2000i i for 
an age of r = 537 Myr and Z = 0.008 ([Fe/H] = -0.38) tfiat 
provides the best match to the observations of Capella (Af, Tj-ff, 
L, and R for both stars). Two different evolutionary states for the 
primary provide similarly good fits: one on the ascending branch, 
and the other on the descending branch prior to the beginning of of 
the core helium-burning phase. The large filled circles and shaded 
boxes correspond to the measured quantities for Capella and their 
uncertainties, while small dots on the isochrone correspond to the 
locations of the best fits. See Table [TOl 

most previous investigators. Although the detailed sur- 
face abundances are not include d in the published ta bles 
for these models, recent work bv lBertelli et al.l ()2008f ) ex- 
tending the same series of calculations does provide typ- 
ical values for some elements. For a 2.5 Mq star similar 
to the Capella primary and a metallicity oi Z = 0.008, 
the predicted ^^C/^^C after the first dre dge-up is 19.0. 
This i s lower than the estimate of 27±4 bv lTomkin et al.l 
(|1976f ). On the other hand, the C/N ratios show much 
closer agreement with theory. The predicted value for the 
primary at the end of first dredge-up is 0.67, which is only 
slightly higher than the measured value of 0.57 ± 0.06. 
For the secondary, which has not yet experienced first 
dredge-up and may therefore be expected to still have its 
initial main-sequence value of 3.27, the measured value 
is in fact very close: 3.30 ± 0.16. 

In order to provide more flexibility in exploring possi- 
ble evolutionary stages for the primary, and also to gain 
better access to the detailed surface com position we h ave 
considered a different set of models bv IClare"3 ()2004[ ) in 
which we can more easily vary not only the overall metal- 
licity but also the mixing length and overshooting param- 
eters. Evolutionary tracks were computed for the exact 
masses we measure. The best match is found for the 
same metallicity as before {Z = 0.008, corresponding 
to [Fe/H] = —0. 33 in these models, which adopt solar 
abundances from iGrevesse fc Sauval Il998l ) and convec- 
tive parameters auL = 1-63 and Oov = 0.20. As indi- 
cated in Figure fT2k . the primary is on the descending 
branch, prior to the core helium-burning phase, and the 
age determined for the two components (553 Myr) is vir- 
tually identical to within 0.2%. Details of the differences 
in Tcff , L, and R for each star are given in Table 1171 
(case A), along with the goodness of fit (x^)- Lithium 
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Fig. 12. — Evolutionary tracks f or Capella A (solid) and B 
(dashed) from the series of models bv lClaretl l |2004l ). for the exact 
masses determined here. The measured temperature and luminos- 
ity of each component are shown with the filled circles and error 
boxes. A thin line connects the measured positions with the best- 
fit location on each track. Four different scenarios are considered, 
focusing on the evolutionary status of the primary star. The pa- 
rameters of each case and the location of the primary are as labeled. 
The asterisk on the primary tracks indicates the position where the 
models give a ^^C/^^C carbon isotope ratio exactly matching the 
measured value. 
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Fig. 13. — Carbon isotope ratio evolution for the primary of 
Capella. (a) ^^C/^^C ratio as a function of age according to the 
Claret model in Figure [I2}i {Z = 0.008, ciml = 1-63, Oov = 0.20; 
case A). The measurement is indicated with the dashed line, and 
the shaded area represents the la uncertainty, (b) Location of the 
bottom of the convective zone in terms of the radius of the star, 
showing the rapid deepening during the first dredge-up. (c) Close- 
up of the ev olutionary track for the measured mass of Capella A 
from Fig. 112b showing the change in the carbon isotope ratio during 
the first ascent. The measured temperature and luminosity along 
with the error box are indicated. An asterisk marks the location 
on the track where the isotope ratio matches the measured value 
of 27 ± 4. 

abundance calculations from these models are not up to 
date and have been found to give surface values that are 
much lower than other calculations. They are not con- 
sidered further, although this has no impact on other 
calculations since Li is a low abundance element with 
no energetic importance. As in the previous models, the 
predicted ^^C/^^C isotope ratio for the primary is found 
to be lower than measured (19.18 versus 27 ± 4, nearly 
a 2(7 difference) . At face value the measured ratio seems 
inconsistent with any position for the primary of Capella 
other than on the first ascent. This is illustrated in Fig- 
ure [T31 The decrease in the isotope ratio shown in the 
top left panel is very rapid and follows closely the deep- 
ening of the convective envelope displayed on the lower 
left (first dredge-up). Values of ^^C/^^C are indicated 
on the evolutionary track in the right panel. The ratio 
changes from 85 to 20 in less than 1 Myr, so there is 
very little leeway to accommodate the measured value, 
which is predicted to occur in this model at a consider- 
ably lower luminosity and somewhat hotter temperature 
than observed, by 2.9a and l.Ocr, respectively (asterisk 
in Figure \Wb ) . The alternative is an error in the mea- 
surement. We discuss this further below. On the other 
hand, the C/N ratios expected from these models are 
quite consistent with the observed values, as seen in Ta- 
ble [HI 

Given the disagreement of the previous models with 
12q^13q^ we explored other possible matches that place 
the primary on the ascending branch, by varying both 
the mixing length and overshooting parameters. One 
such solution has Z = 0.004 (corresponding to [Fe/H] 
= —0.63), auL — 1-44, and a reduced overshooting of 
Q^ov = 0.10. The ages inferred for the two components 
agree to within 1.0%, and have a mean of 464 Myr (see 
Figure [l2b) . The carbon isotope ratio for the primary 
is underpredicted by about the same amount as before 
(Table [TTl case B), and the location on the track that 



matches the measured ^^C/^'^C value (asterisk) is only 
marginally closer to the measured luminosity and tem- 
perature of Capella A (2.7ct and O.Bcr discrepancies, re- 
spectively) . A similar fit with the primary also on the as- 
cending branch can be obtained with qml = 1-30 and no 
overshooting (aov = 0.00), but at an even lower metallic- 
ity of Z = 0.002 that is probably unrealistic. The mean 
age in this case is 400 Myr, and the individual ages differ 
again by only 1.0% (Table [HI case C, and Figure [Efc) . 
The luminosity and temperature at which the primary 
model predicts a carbon ratio matching the observation 
are 2. Oct and 0.9ct away from the measured values, re- 
spectively. Setting aside the ^^C/^^C discrepancy for the 
moment, it is interesting to note that the measured C/N 
ratio for the primary appears to rule out all scenarios 
with that star on the ascending branch. While the pre- 
dictions for (C/N)b are still consistent with the obser- 
vations, the expected (C/N) a values are a factor of ~3 
too high (Table flTl cases B and C). Thus, the chemical 
evidence for Capella seems contradictory regarding the 
viability of the first ascent scenario. 

We investigated also possible fits to the Claret mod- 
els with the primary in the core helium-burning phase. 
While previous researchers have had no difficulty in find- 
ing satisfactory fits, in this case we are unable to achieve 
a match to M, R, Tcff, and L as good as in the scenarios 
discussed above. Evolutionary tracks with Z — 0.008, 
OiML — 1.40, and ctov = 0.20 yield ages for the stars 
that differ by 16.4% (mean age of 589 Myr). The carbon 
isotope ratio for the primary is not reproduced (see Ta- 
ble [HI case D, and Figure [T^). but the C/N fractions 
are consistent with the measurements since the primary 
is well past the first dredge-up. 

None of the stellar evolution models considered so far 
include the effects of rotation. Although the current ro- 
tation rates for Capella A and B arc slow or moderate, 
this was not the case for the main-sequence progenitors, 
which were late B- or early A-type stars spinning at typ- 
ical equatorial velocities of 150-200 kms~^. Evolution- 
ary tracks with rotation were computed for both com- 
ponents following Clare3 (|1999l ) for the same set of pa- 
rameters considered in case A above, under the assump- 
tion of rigid-body rotation. The initial angular rotation 
rate Wrot was set by requiring an approximate match be- 
tween the predicted equatorial velocities at the current 
age and the measured values (de-projected from vsini 
by assuming the spin axes are aligned with the orbital 
axis; see below). However, given that the rotation of the 
primary appears to be synchronized with the orbital mo- 
tion (S I9.2|) . one may expect tidal forces to have had some 
effect on that star whereas the secondary is farther from 
synchronization, and more likely to have evolved as an 
isolated star. Thus we only required a match to the rota- 
tion of the secondary, and we adopted the same value of 
Wrot for the primary on the grounds that the masses are 
the same within 1%, so all properties including rotation 
must have been very similar to begin with. The differ- 
ences compared to non-rotating models are most obvious 
on the main sequence, where rotation renders the tracks 
some 500 K cooler, but are smaller toward the later evo- 
lutionary stages. We find that the best match between 
the rotating models and the observations is obtained for 
parameters nearly identical to those found earlier, at an 
age only 5 Myr older. The details of this comparison are 
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Fig. 14. — Evolution of the rotation rates of Capella A (solid 
line) and B (dashed line) relative to their values on the ZAMS, 
according to Claret models with Z = 0.008, oml = 1.63, and Oqv = 
0.20 (case A). The inset shows the ratio of these curves (secondary 
divided by primary) along with the measured ratio of the projected 
rotational velocities (i)b sini/v\ sini = 7.1 ±1.5; shaded horizontal 
band). The theoretical ratio agrees with the measurement at the 
evolutionary age of the system (shaded vertical band spanning the 
primary and secondary ages from Table I17I I . This suggests tidal 
forces within the binary, which are not accounted for here, are 
weak. 

included in Table flTl case E. The older age produced by 
these models compared to case A is reminiscent of the 
effect of overshooting, which extends the main-sequence 
lifetime. We also find that, although we did not require 
a match, the predicted rotational velocity of the primary 
reproduces the measured value exactly. This is an in- 
dication that the effects of tidal forces in damping the 
rotation of the star are relatively weak, given the long 
orbital period and relatively small relative radii of the 
components {Ra/ci ~ 0.075 and Rb/o- ~ 0.055, respec- 
tively) . Another way of looking at this is to examine the 
evolution of the rotation rate for each star due solely to 
changes in the moment of inertia with time, under the as- 
sumption of rigid-body rotation and conservation of an- 
gular momentum without including tidal forces. These 
rotation rates relative to their initial values on the zero- 
age main sequence (ZAMS) are depicted in Figure [HI for 
both components as a function of time, for the models 
with ttML — 1.63. The inset shows the ratio of these 
two curves near the present age of the binary, which 
should correspond to the actual ratio of the rotational 
velocities if the stars started out with the same value on 
the ZAMS and if tidal forces are unimportant. There 
is in fact good agreement with the measured ratio of 
ub sin i/wA sin i = 7.1 ± 1.5 (horizontal shaded area) at 
the evolutionary age determined with these models (ver- 
tical shaded area). 

The compari sons made above aga inst the stellar evo lu- 
tion models bv lGirardi eFan (|2000f ) and lClaretl (l2004[ ) do 
not constitute very strong tests of the predictive power 
of theory because a number of parameters have been 
adjusted to produce a good fit to the observations of 



Capella, namely, the overall composition, the mixing 
length, and overshooting. We note that in both mod- 
els the metallicity that seems to be favored {Z = 0.008) 
is actually in very good agreement wi th the spectroscopi c 
value (§[7l) derived from the work of McWilham' (l990), 
despite our concerns about possible biases in those mea- 
surements. Even if we accepted the latter as accurate, 
and had imposed that value on the models thereby elim- 
inating Z as a free parameter, oml and Oov are still ad- 
justable quantities. This is unavoidable because of the 
phenomenological way in which convection and mixing 
are treated in these "standard" models. 

We have thus considered a third se t of models based on 
the TYCHO stellar ev olution code (|Young et al.ll2001al : 
lYoung fc Arnetti l2005f ) that implements somewhat dif- 
ferent prescriptions for some of the physical processes. 
Predicted surface abundances are available for all the el- 
ements in the network, including lithium, and rotation 
has not been considered for the present work. The mod- 
els calculated with the most recent version of this code, 
TYCHO 8, have a partial implementation of results from 
a new, non-l ocal hydrodynamic description of convection 
presented in lArnett et al.l ()2009f ). The key results imple- 
mented in this version are the boundary conditions for 
convective zones. Convective stability is evaluated based 
on the Richardson number {Ri = N'^ / {du / dr)'^) . This 
compares the potential energy in stratification as mea- 
sured by the Briint-Vaisala frequency TV versus the ki- 
netic energy in shear (including turbulent velocities). Us- 
ing the Richardson criterion results in larger convective 
zones since thermodynamically stable regions in mixing 
length theory can be hydrodynamically u nstable. The 
compl ete velocity integration described in lArnett et al.l 
(|2009( l has not yet been implemented for this work. 
Convective velocities are calculated according to mix- 
ing length theory with a mixing length equ al to the 
depth of the convection zone, according to lArnett et al.l 
((20091). This eliminates Oml as a free parameter. For 
very deep convective zones the simulations suggest there 
may be a saturation in the mixing length at a value of 
roughly 1.5-2.5 as dissipative processes become impor- 
tant. We place this saturation near the lower limit at 
1.64. The models also include wave-generated mixing 
from i ntermittent turbul e nce an d slow circulation accord- 
ing to lYoung fc Arnettj (|2005D . For most of the star's 
history the convection is relatively shallow {Iqz < 2/ip, 
where Iqz is the depth of the convection zone and hp 
is the pressure scale height) except in the photosphere. 
In this regime the convective description has no free pa- 
rameters. On the RGB, however, the convection zone 
has a very large density gradient over its depth. This 
has two effects. First, significant radiative losses from a 
fluid parcel can occur over a single convective traverse. 
Secondly, the sound speed drops with increasing radius 
and decreasing density. Initially strongly subsonic flows 
can generate shocks in the outer regions of the convec- 
tive zone causing more dissipation. Unfortunately, we do 
not yet have a satisfactory quantitative theoretical de- 
scription of this process, so we are forced to parametrize 
it. We introduce a geometric parameter ^ml (compare 
to the geometric parameter introduced by Bohm-Vitens^ 
(IT958h for mixing length) which describes the aspect ra- 
tio of a convective fluid parcel. The value of the Bohm- 
Vitense parameter gives a surface-to-volume ratio less 
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Fig. 15. — Evolutionary tracks for the me as ured m asses 
of Capella from the models of lYoung et al.l I l2001al ) and 
lYoung fc ArnettI | |2005D . Two different scenarios are shown, fo- 
cusing on the evolutionary state of the primary. 

than that of a sphere. If we compare our parameter to the 
mixing length context, values of 1 to 1000 are a reason- 
able range to account for true geometric differences for 
large surface area turbulent flows, radiative losses over 
large temperature gradients, and hydrodynamic dissipa- 
tive processes. 

Using these models in the same way as before, the 
best match to the observations for Capella was found 
for evolutionary tracks with an overall metallicity of 
Z = 0.4:5Zq (in whi ch th e solar composition is adopted 
from lXsplund et al.ll2005( ). corresponding approximately 
to [Fe/H] = —0.35. This is very similar to the composi- 
tion that gave the best fits in the previous sets of models. 
For the secondary we used guL = 1 since it is only begin- 
ning to develop a deep convective envelope, and we are 
using a parametrization instead of a quantitative physical 
description. A modest increase may be supportable and 
provide a better fit. For the primary a value of 5ml — 100 
provides a reasonable match to the global properties M, 
R, Teff, and L with the star on the red clump fTable [T51 
case A), giving a mean age for the binary of ^593 Myr 
and an age difference between the components of about 
6%. A value for gML of 1 gives a good fit for the primary 



on the first descent of the RGB, but a poorer age match 
with the secondary (see Table [T51 case B). An interme- 
diate value of guL = 10 yields a less satisfactory agree- 
ment, also on the first descent. The above comparisons 
are shown graphically in Figure [T51 The C/N abundance 
ratios predicted for the secondary in these two scenarios 
are quite consistent with the measured value, while those 
for the primary are more discrepant (see Table [18]). 

The ^^C/^'^C ratios expected for the primary from the 
TYCHO models (19.0-19.6, depending on the evolution- 
ary stage) are still significantly lower than the measured 
value of 27 ± 4, as found b efore with the Clare t models 
(18.9-19.2) as weU as the iBertelh et all ()2008D models 
(~19 .0). Other calcula tions indicate very similar values 
(e.g.. ICharbonn el 19941. It would appear, therefore, that 
current theory systematically underestimates the carbon 
isotope ratio in this case, regardless of where the primary 
is located on the H-R diagram. Unfortunately there is 
only a single measur ement of ^^C/^^C for Capella A, by 
iTomkin et all (fl976l ). A comparison of the results from 
that work for other stars that have also been measured by 
others seems to indicate some systematic differences, but 
it is difficult to te ll which values are co rrect. For example, 
measurements bv lTomkin et al.l (|1976l ) of the four Hyades 
giants 7 Tau, S Tau, e Tau, and 0^ Tau, which have 
masses only slightly smaller than Capella, give -'^^C/^fC 
ratios of 19, 23, 22, and 20, respectively, while iGilrovl 
(|1989D reported 24, 26, 26, and 27, aU systematically 
higher. It is possible that errors in the determination 
of the microturbulent velocity £t can explain these dif- 
ferences, since Keller et al. (2001) and many others have 
shown that the ^^C abundances are quite sensitive to 
because the ^^CN lines from which they are often mea- 
sured are typically close to saturation. A cursory exami- 
nation of measurements for other giant stars by different 
authors indicates similar discrepancies, which leads us 
to believe the uncertainties in some of these studies may 
be underestimated. A new measurement of the ^^C/^'^C 
ratio for Capella A would be extremely helpful to settle 
this issue. 

TablefTSIshows that the surface lithium abundance pre- 
dicted for the secondary of Capella from the TYCHO 
models is in agreement with the measurement, within 
the uncertainties, and is not much reduced from its main- 
sequence value. The primary lithium abundance, on the 
other hand, is larger than observed by about 3a, re- 
gardless of the value of guL- According to these models 
log7V(Li) is not expected to dip much below 1.9 even af- 
ter the phase of core helium burning, as the star evolves 
up the asymptotic giant branch. Thus, it appears lithium 
burning in Capella A has been more efficient than pre- 
dicted by theory during the first dredge-up, or there 
was additional burning afterwards. Unlike the case of 
^^C/^'^C, there are several independent measurements of 
the Li I A6708 line for the primary, and they display 
rather good agreement. A possible explanation for the 
discrepancy is that in the current TYCHO models rota- 
tional and convective/wave sources of shear at the con- 
vective boundary are not correctly coupled. Mixing at 
the base of the envelope in the models will be less than 
that of the real star, and this could result in discrepancies 
in the abundances of low-temperature burning products. 

As seen from the comparisons above with three differ- 
ent sets of evolutionary models, we are unable to find 
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an evolutionary state for the binary that gives complete 
consistency with all absolute dimensions for both stars 
as well as all chemical indicators, simultaneously. Par- 
tial matches are possible, indicating either deficiencies in 
our understanding of stellar evolution, observational er- 
rors in some quantities, or both. Other evidence of the 
evolutionary status of Capella has occasionally been con- 
sidered by some authors, including activity indicators in 
the UV and X r ays. For example, b ased on similarities 
with other stars, lAvres et al.l ()1983[ ) estimated that the 
primary may be responsible for as much as half of the 
total 0.1-0.4 keV soft X-ray coronal emission from the 
system, even though it only accounts for a small frac- 
tion of the far-UV emission characteristic of the chromo- 
spheric layers. A comparison with the generally weaker 
coronal X-ray levels of clump giants led them to con- 
clude that Capella A is crossing the Hcrtzprung gap for 
the first time, and is yet to go through the first ascent. 
Although this significant level of X-ra y emission was sub- 
sequently confirmed by others (e.g.. iLinskv et ahl 119981 : 
lYoung et al . 2001b), this picture seems inconsistent with 
some of the chemical indicators discussed above. Ad- 
ditionally, clump giants are known to display a wide 
range of coronal emission levels (see, e.g.. lPizzolato et"al] 
12000) , perhaps indicative of activity cycles analogous to 
that of the Sun.^^ Virtually all UV and X-ray studies 
after the mid 1980s have implicitly assumed the primary 
of Capella is a core-helium burning star, on the grounds 
that it would be very unlikely to find both components 
of a binary in a very rapid state of evolution such as the 
Hcrtzprung gap or the first ascent, given that the sec- 
ondary is already known to be a first crosser. This could 
only happen if the mass ratio were very near unit y, and 
that d id not appear to be the case until now.^^ lAvresI 
(|1988f ) proposed the primary may be in the early stages 
of core-helium burning, while other giants that are even 
more active such as /3 Cet and the Hyades stars 7 Tau 
and 9^ Tau were considered to be in a more advanced 
clump phase ("blue loop"). Qualitative assessments of 
the evolutionary state based on the overall activity level 
remain problematic, and are not easily testable against 
models. 

We conclude this section by pointing out that the de- 
tached nature of Capella, and therefore its value for test- 
ing stellar evolution theory, is beyond doubt. The size 
predicted for the primary at the tip of the giant branch 
from any of the models considered above is less than 
35 i?0, and would still fit comfortably inside its criti- 
cal Roche surface (~60 i?©). Thus, even if the primary 
were already in the clump, mass transfer by Roche-lobe 
overflow would have been avoided. 



9.2. Comparison with tidal theory 

There may in fact be evidence that C apella itself is undergoing 
such cycles. Kat sova fc Scherbakovl l|1998l ) detected long-term vari- 
ations in the equivalent width of the chromospheric He I 10830 A 
line assoc iated with the primary, wit h a possible period of roughly 
11 years. IRaassen &: Kaastral I I2007I) found a similar variability in 
X rays in the flux of the coronal Fe XXI 128.73 A line, also asso- 
ciated mainly with the primary, coincidentally with an apparent 
period also near 11 years. 

The most recent determination, bv I Barlow et al.l ()1993l l. gave 
a mass difference of 5.2%. In our analysis it is less than 1%. 



Capella is a particularly favorable case to study the 
rotational state of the components. In addition to the 
measures of wsini and the radii, the rotational peri- 
ods Prot of both stars have been measured as well. 
IShcherbakov et aL (1990) found that the radial velocity 
of the chromospheric He I A10830 absorption line follows 
the velocity of the photospheric lines of the primary, un- 
ambiguously identifying this activity feature with that 
star. The equivalent width of the line was also seen to 
vary with a period indistinguishable from the orbital pe- 
riod, which wa s interpreted as the rotationa l signature 
of the primary. iKatsova fc Scherbakovl (jl998f ) confirmed 
this with additional data, and measured a rotation period 
of 103.97 days. They were also able to detect a weaker 
component of the He I line associated with the secondary 
star. The radial velocity of this feature again varied in 
phase with the photospheric velocities for Capella B, but 
its equivalent width was seen to vary with a period of 
8.25 days, repres e nting the rotation of the secondary. 
IStrassmeier et all (|2001t ) detected the rotational signa- 
ture of both stars in their Ha photometry, with periods 
of 106±3 days for the primary and 8.64±0.09 days for the 
secondary. These observations make it very likely that 
the primary's rotation is effectively synchronized with 
the orbital motion, while the secondary is spinning more 
than 12 times more rapidly. For the purpose of this work 
we adopt a compromise value of Prot = 8.5 ± 0.2 days for 
Capella B. 

With these periods and the measured radii, the pre- 
dicted rotational velocities projected along the line of 
sight are 3.92±0.19 kms^^ and 35.4± 1.5 kms^^ for the 
primary and secondary, respectively, under the assump- 
tion that the spin axes are aligned with the orbital axis. 
These predictions may be compared directly against the 
measured wsini values in Table [HI We note, however, 
that most of those determinations — including our own 
— have ignored the additional broadening that comes 
from macroturbulence. This can be significant for gi- 
ants and in our case is particularly important for the 
primary, which is a slow rotator. The only two determi- 
nations of wa sin i in the table tha t have accounted for 
macr oturbulence (|Fekel et al.lll986HStrassmeier fc Fekel 
Il990h are both smaller than ours, and average 4 kms . 
A giant star with the spectral type of Capella A is ex- 
pected to have a macroturbulent velocity Crt between 5 
and 6 kms~^ (see, e.g., [Gray 1982), whereas the syn- 
thetic templates we used in S I2.1I to determine v sin i 
have Crt = 1 .5 kms~^. A cor rectio n for this differ- 
ence f ollowing I Fekel et ahl (|1997D and iMassarotti et al.l 
(|200l brings our v sin i value for the primary down from 
6.5 kms~^ to about 5.0 knis^^. For the secondary this 
effect hardly matters, but a similar correction assuming 
~ 7 kms~^ leads to WBsini — 35.6 kms~^. These 
are the values we adopt in the following (Table [T5|) . The 
agreement with the predicted velocities, which is well 
within the errors, can be considered as an indication that 
the absolute radii based on angular diameters are accu- 
rate. 

The binary system of Capella resembles in many re- 
spects that of the 76-d ay period eclipsing system TZ For 
() Andersen et al.l[l991h . which is also composed of giants. 
Both have a circular orbit, and one component in syn- 
chronous rotation while the other is super-synchronous. 
This puzzling situation is of course explained by tidal 
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forces, which increase exponentially as the more massive 
primaries evolve first and expand in size. The result is 
a sharp drop in their rotation rates, and the damping of 
the orbital eccentricity. We have computed the critical 
times of rotational synchronization and orbital circular- 
ization for Capella according to the turbu lent d i ssipa- 
tion and radiative damping mechanisms by IZahnI (|1992l . 
and references therein), as well a s the hydrodynamical 
mechanism of iTassoul fc Tassoull ()1997l and references 
therein). Note that we calculate actual times rather than 
timescates, since the latter would be completely meaning- 
less for giant stars such as Capella that have altered their 
structure drastically since the ZAMS, changing from hav- 
ing radiative to convective envelopes. Thus we have in- 
tegrated the differential equations for the changes in the 
orbital eccentricity and rotation rates alo ng the evolu- 
tionary trac ks following the prescriptions of lClaret et al.l 
(|1995l [T997[) . to properly account for the stellar proper- 
ties during the radiative and convective phases of evolu- 
tion. The results obtained using the Claret non-rotating 
models for Capella A and B for Z = 0.008, oml = 1-63, 
and Q!ov = 0.20 (case A) are illustrated in Figure [TBI 
In this diagram the radii are plotted as a function of 
age, and the times of synchronization and circulariza- 
tion are indicated along with the evolutionary age of 
the sy stem (f = 5 53.2 Myr). For the primary the the- 
ory of IZahnI ()1992f ) gives Tsync,A — 545.8 Myr, and for 
the secondary Tsync,B — 559.8 Myr, both of which agree 
with the observation that the primary already appears 
to have its rotation synchronized with the orbital mo- 
tion while the secondary is clearly asynchronous. Ac- 
cording to these models the spin of each star becomes 
locked with the orbital motion when the components 
reach their maximum size at the tip of the giant branch 
(Figure [TBI top). Circularization occurs somewhat later, 
at Tcirc — 563.6 Myr. Strictly speaking, this is larger 
than the age of the binary and the prediction is thus for- 
mally inconsistent with the apparent circularity of the 
orbit, although perhaps not by much considering the ap- 
proximations in tidal theory^^. Similar calculations for 
the mechanism of Tassoul & Tassoul (1997) are shown in 
the bottom panel of Figure I16i and indicate formal con- 
sistency with the observations for all three critical times. 

We have repeated the calculations according to both 
tidal theories for the other scenarios considered earlier 
based on the Claret models, and the agreement with 
the observations varies. The results are summarized in 
Table [T^ With one exception ( case A and the mech- 
anism of ITassoul fc Tassoull [1997[ ). in no other scenario 
are all three critical times strictly consistent with real- 
ity, but the predictions are so close to the best-fit evo- 
lutionary ages that the results are very sensitive to the 
details of the models and the calculations. For example, 
the mere addition of rotation to the first set of evolu- 
tionary tracks examined above is enough to change the 
conclusions entirely regarding tidal theory, even though 

^"^ In particular, the tidal equations used here are valid only 
for small departures from synchronization, and do not include the 
effect that magnetic fields in these active stars could have on an- 
gular momentum transport. The secondary of Capella is rotating 
some 12 times faster than synchronous, so the predictions should 
be taken with caution. An additional approximation in our calcula- 
tions is that the concomitant changes in the orbital semimajor axis 
are ignored, although this is likely to be a small effect for Capella. 
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and circularization (Tdrc) for the components of Capella from tidal 
theory, shown in a diagram of radius versus age. Solid and dashed 
lines represent the predicted radii for the primary and secondary, 
respectively, according to Claret models with Z = 0.008, Qml = 
1.63, and «ov = 0.20 (case A). The measured radii are shown as 
horizontal shaded regions, and the evolutionary age of the system 
is indicated with a vertical shad ed area. Calculations according to 
the tidal theory of Zahn' l'1 992l. and refere nces therein) are shown 
at the top, and [Tassoul & Tassoull l|1997l . and references therein) 
at the bottom. 



it hardly changes the stellar properties for the best fit. 
Thus, it is not possible to establish which tidal theory is 
more appropriate for Capella. We note, in this connec- 
tion, that the theory of Tassoul has often been found to 
be too efficient, at least for main-sequence stars, giving 
critical times that are occasionally 1-2 orders of magni- 
tude shorter than those by Zahn. In Capella that is not 
the case. Given that the critical times from these two 
theories tend toward each other in the limit of long pe- 
riods, this again is an indication that tidal forces in this 
binary are relatively weak. 

Another consequence of tidal forces in binaries is that 
the spin axes of the stars tend towards alignment with 
the orbital axis. Because the timescale for this is typi- 
cally much shorter than for the other two effects, at least 
for main-sequence stars, it is almost universally assumed 
that the axes in binary systems are aligned. The sim- 
plifying advantage is that one can then set the inclina- 
tion of the rotation axes equal to the orbital inclination, 
as we ourselves have done above, for example to infer 
equatorial rotational velocities. We note however, that 
this assu mption may not always hold (see, e.g., the case 
of A Vir: IZhao et al.|[2007f ). The timescale for alignment 
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depends strongly on the ratio between the orbital and ro- 
tational angular momenta. This ratio is large for Capella 
because of the long orbital period, so the timescales for 
alignment are close to those for synchronization. As an 
example of such calculations, we use the rotating Claret 
models with Z = 0.008, q;ml = 1-63, and aov = 0.20 
(case E), and find t hat th e alignment times according to 
iTassoul fc Tassoull ()199/t ) are 555.9 Myr for the primary 
and 563.6 Myr for the secondary. These are again so near 
the evolutionary age of 558 Myr for this scenario that it 
is difficult to draw general conclusions. At face value the 
primary is formally predicted to have reached alignment, 
while the secondary is still a few Myr fro m reaching it . 
The critical times for alignment according to lZahnI ()1992f ) 
are 562.3 Myr and 570.2 Myr, both of which are formally 
greater than the age, implying that the spin axes are not 
necessarily aligned. In Capella there is an observational 
constraint on this prediction, given by the rotation pe- 
riods, wsini values, and absolute radii determined for 
the two components. As indicated earlier in this section, 
there is very good agreement (particularly for the more 
discriminating secondary) between the measured v sin i 
values and the predicted rotational velocities, when the 
latter are projected along the line of sight using the well- 
known orbital inclination angle. This is evidence that 
the orbital axes are in fact closely aligned with the orbit, 
and we conclude that tidal theory does not clearly agree 
with the observations in this instance. 

10. CONCLUDING REMARKS 

With our new spectroscopic observations of Capella 
and a global orbital solution that uses all available as- 
trometric and radial-velocity measurements, the preci- 
sion of the most important of the physical properties of 
the stars — their masses — has been improved threefold 
compared to previous values, and now stands at 0.7% for 
the primary and 0.5% for the secondary. These are at the 
level of the best available determ inations from do uble- 
hned eclipsing binaries (see, e.g.. lAndersenI [T99l[ ). As 
others have pointed out, the mass ratio is perhaps more 
important in this case since it sets constraints on the 
differential evolution of the stars; that quantity is now 
known to even higher precision (0.36%). Based on the 
tests and crosschecks documented in this paper, it is our 
belief that the accuracy has also improved significantly. 
Our dynamical masses differ from earlier results by 9% 
and 5%, respectively, which has an enormous impact on 
the evolutionary state of giants such as these. 

Our goal here has been to bring to bear all relevant 
information on the physical properties of Capella for a 
comprehensive discussion of the system in terms of its 
evolution. In addition to the masses, we report signifi- 
cantly improved determinations of the absolute radii, the 
effective temperatures, the luminosities, the rotational 
state of the stars, and the chemical composition, in- 
cluding the overall [m/H] index, the lithium abundance, 
and the ^^C/^^C and C/N ratios. We know of no other 
evolved system with so much high-quality information 
available. Sanity checks have been performed, whenever 
possible, to assess the accuracy of these quantities. Our 
expectation at the outset was then that we should be 
able to pinpoint the evolutionary state of these giants, 
particularly the primary, and validate current models of 
stellar evolution that seem to work so well on the main 



sequence. Instead, we find that there are now so many 
constraints that no single model is able to satisfy them 
all simultaneously, even when adjusting some of the free 
parameters in theory. 

Three types of stellar evolution models have been 
consid ered, two based on s tandard cony ection prescrip- 
tions (jOirardi et all 120001 : IClaretl |2004[ ) and one with 
an alternative form ulation (TYCHO evolution code; 
lYoung et al.ll200l"al ). We have also explored the role of 
rotation. In most cases it is possible to obtain reason- 
ably good matches to M , R, L, and T^s for both compo- 
nents at a variety of evolutionary states for the primary, 
of which the clump does not particularly provide one of 
the best matches. This is in contrast with most previ- 
ous studies, which overwhelmingly place the primary in 
the core-helium burning phase, and is largely due to the 
fact that the mass ratio is now much closer to unity than 
previously believed (within 1%, versus 5% previously). 
This allows the "evolutionary gap" between the com- 
ponents to narrow considerably. These four properties 
alone are therefore insufficient to discriminate between 
different models or different states for Capella A. 

Chemical constraints are seen to be particularly pow- 
erful in this case. We find that the overall metallicity 
[m/H] required for the models to match the properties of 
the two stars is in excellent agreement with the spectro- 
scopic determination. However, the predicted ^^C/^'^C 
values for the primary are systematically too low (at the 
~2cr confidence level), if the measurement is to be taken 
at face value, and would require that component to be 
near the base of the RGB. This is inconsistent with the 
well-determined temperature and luminosity (or radius). 
Errors in the rather delicate measurement of the carbon 
isotope ratio cannot be ruled out, and a new determina- 
tion would be highly beneficial. Similarly, models tend to 
overestimate the apparently well-measured lithium abun- 
dance of the primary (discrepancies at the ~3f7 level) , al- 
though they do match that of the secondary, which is es- 
sentially the initial value. The observed C/N ratios only 
permit the primary star to be on the descending branch 
prior to the helium-burning phase, or at the clump. As- 
cending branch scenarios are ruled out at the many-cr 
level. A location on the descending branch would in fact 
be more likely than on the first ascent, given that the 
descent takes about 3 times longer (see footnote [T^ . 

It may well be that the primary is indeed a clump gi- 
ant, but the fact remains that current stellar evolution 
models fail to satisfy all of the observational constraints 
when considered jointly, for both components at the same 
time, and at a single age, whether the primary is in the 
core-helium burning phase or at any other location. Pre- 
dictions are still very strongly dependent on the physics 
in the models, underscoring our incomplete understand- 
ing of many of these processes. The ^^C/^^C measure- 
ment suggests that convective mixing in the primary is 
less efficient than predicted, while the lithium determi- 
nation appears to shows the opposite. It is also quite 
possible that the masses, as precise as they already are, 
need to be improved further by an order of magnitude 
given the sensitivity of the predictions. If that is the case, 
then both better spectroscopy and better astrometry will 
be required, since the latter contributes a non-negligible 
^--^0.3% to the present mass errors. 

Tidal theory does not seem to fare much better. The 
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circular orbit of Capella along with the evidence for syn- 
chronous rotation of the primary and asynchronous rota- 
tion of the secondary pose strong constrain ts that are in 
gen eral not quite met by the mechanisms of lZahnI ()1992[ ) 
and'Tasso ul fc Tassoull |1997). Similarly with the align- 
ment of the rotational axes with the orbital axis, which 
is observationally constrained by the measurements of R, 
vsini, and Prot- To be fair, however, the critical ages for 
circularization, synchronization, and alignment inferred 
from these models for the various scenarios considered 
are typically close to the evolutionary age in each case, 
so the results may be sensitive to the details of the cal- 
culations and are not conclusive. 

Despite these challenges, much has been learned about 
Capella in its more than 100 years of observational his- 
tory. It has been a favorite target of virtually every space 
facility and has provided, among many other results, very 
important insight into the structure and other properties 
of the coronae in active stars. Because it is so bright, 
nearby, and easy to observe, Capella remains a very im- 
portant benchmark for testing our understanding of stel- 



lar physics in the advanced stages of evolution. 
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APPENDIX 

THE GREENWICH VISUAL OBSERVATIONS OF CAPELLA 

The relatively long orbital period of Capella from the spectroscopic studies of iNewalj (|1899[ ) and ICampbelll ()1901[ ) 
immediately suggested to investigators from the beginning of the 20th century the possibility that it might be detected 
as a visual binary, with a separatio n somewhat u nder O'.'l, according to the best estimates at the time. Foreseeing 
the potential of interferometry. Miller Barr (1900) suggested it should be resolvable with Michclson's "interference 
apparatus" on a large telescope, although this would only come to happen two decades later. In the meantime, visual 
binary observers lost no time in attempting the measuremen t, even though predictions indicated the separation would 
be at or below the resolution limit of available telescopes. lAitkenI (|1900r ) described five separate attempts with the 
Lick Observatory 36-inch Clark refractor (resolving power ~0'.'15), and on only one of those dates (1900 February 24) 
did he indicate that once or twice a slight elongation of the image was suspected, but that repeated measurements of 
the position a ngle showed a spread of some 60°, which he considered large enough to be a sign that no real elongation 
was observed. iHussevI () 19001 . [19011 ) used the same telescope to look at Capella on multiple occasions, but also reported 
the images to be essentially round. 

At about the same time, visual binary observers at the Greenwich Observatory took an interest as well, and reported 
seeing the image elongated with the 28-inch refractor (resolving power ^0'.'19) on 1900 April 4. They even recorded 
a position angle measurement for Capella. Over the next two years they continued observing the star regularly and 
making measurements of 9, and on a few occasions they even ventured to estimate p. More than 100 p ositiq r i-angl e 
measurements were made by 11 different observers, including the As tronomer Royal (jChristiel Il900al lbl. 1190 21. 119031). 
These observations have been viewed with suspicion ever since (e.g., iBurnhanil 1 1 906f ) . As pointed out by [Ashbroo^ 
(|1976f ). it is difficult to see how "a smaller instrument at sea level could consistently do what the 36-inch Clark refractor 
on Mount Hamilton found impossible". No ot her visual observati ons elsewhere appear to have been successful; the 
only other recorded attempts, decades later bv I WilsonI (|1939l [l94lh . saw no elonga tion in Capella. 

The most remarkable fact about the Greenwi c h obs ervations, as discussed by lAshbrooa (|1976( ). is their overall 
agreement with the astrometric orbit of iMerrilH (|1922| ). which was not published until twenty years later. Merrill 
himself had compared the Greenwich measureme nts w ith his orbit, but other than noting some relatively minor 
discrepancies, he made no judgement about them. iGoosI (fT90S) had also discussed these measurements in connection 
with his own spectroscopic observations, but this was also long before the astrometric orbit was known. Aside from 
these two studies, the Greenwich observations appear to have been largely ignored. 

The refined orbital solution derived in the present work enables us to examine the accuracy of these measurements 
more closely. The circumstances of each observation were reconstructed based on published information.^^ Most of 
the position angle estimates were made during the daytime, some at night, and a few during twilight, and all of them 
at hour angles greater than about 2^ hours (east or west) because of pointing limitations. A group of observations 
from 1900 April 23 to May 11 were considered bv I Christie (|1900bf) to be influenced by color effects, and we have not 

^® In such cases the binary nature of the object is often obvious from the elongated images, and experienced visual observers can still 
make rough estimates of the angular separation, especially if the components are of similar brightness, as in Capella (see, e.g.,[Aitkcn 1964). 
Even if estimates of p are not possible, a meaningful measure of the direction of this elongation (position angle) can sometimes be made 
for systems that are considerably closer than the nominal resolving power. In any case, under these circumstances, the result can depend 
significantly on the observer's mental disposition, not to mention the observing conditions. 

One misprint was corrected for the observation on 1900 November 9 {Christie 1900b|). The time of observation was adjusted to 
correspond to the hour angle also reported for the measurement. In all other cases we found reasonable agreement between the times 
(GMT) and hour angles listed, within about 2 hours. 
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Fig. 17. — Greenwich estimates of the position angle of Capella made from the elongated images. The different symbols indicate whether 
the observation was made during the daytime, twilight, or nighttime, (a) 9 as a function of time, together with the computed orbit from 
the solution in this paper (which excludes these measurements); (b) O — C residuals of the d measurements from the computed orbit; (c) 
Airmass of each observation. 

considered them here. 

The first peculiar coincidence is that the observers correctly recorded the posi tion angle of Capella to be decreasing 
with time (retrograde motion on the sky), which, as noted by lAshbrookl ()1976[ ). they could not have known in 1900 
since the only information available to them was that provided by the spectroscopic orbit. Secondly, the observations 
correctly placed the ascending node in the first quadrant as opposed to the third. Given the knowledge of the 
spectroscopic orbit, this can only be determined unambiguously from visual o bservations and if a difference in brightness 
is detected, which the Greenwich observers apparently saw, as reported by IChristi3 (jlQOObI ). We now know that this 
difference is quite small in the visible, and as pointed out earlier, observers many decades later were still struggling 
with the difficulty of quadrant determination (see §[S]). Furthermore, the actual position angle of the node implied 
by these observations is very roughly the correct value. The modern determination gives 17 = 40?4, and the estimate 
reported by Christie (1900b) after the first year of observation is = 80°. However, this value may be influenced by 
the fact that the Greenwich observers allowed the orbit to be eccentric, and in fact obtained a small but apparently 
significant value of e = 0.05, whereas we now know the orbit is essentially circular (see §|4]). Repeating the calculation 
for a circular orbit, and holding the period and epoch of nodal passage fixed at their known values, we obtain a revised 
estimate of = 57° from the full sample of original Greenwich measurements, which is even closer to the modern 
value. Finally, the orbital inclination angle from this revised solution is i — 152°, also not far from the correct value 
of 137?2. 

In Figure 117b we show the position angle measurements as a function of time, compared with the prediction from 
our own orbital solution. They cover 7 full cycles of the binary. The second panel displays the O—C residuals. A 
decreasing trend with time is clearly visible, although the magnitude of the residuals is quite modest considering the 
difficulty of the measurements, which can easily have errors of tens of degrees. The bottom panel (airmass versus time) 
shows that during the first few months the observers often made the measurements at large airmasses, whereas later 
on they observed closer to the zenith. It is also during those first few months that the residuals display the largest 
positive trend. 

The agreements described above with the modern orbit of Capella are most intriguing. While it is certainly true 
that these measurements were made pushing the limits of what is possible for the visual technique, and that some 
of the most experienced observers of the time failed repeatedly to see any convinci ng elongat i on in the star's image 
even using larger telescopes, the Greenwich observations are not easily dismissed. lAshbroo"3 ()1976l ) suspected that 



IMerrilll 1)19221 ') had also noticed positive residuals averaging 
might be decreasing with time (see §|3]l. 



-18° from his orbit, and speculated that the position angle of the node 
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some of the measurements could have resuhed — consciously or subconsciously — from foreknowledge of what the 
motion of Capella should be, based on the spectroscopic orbit that had just been published. The observers would 
certainly have known the average daily rate of change of the position angle, but nothing about variations in that rate 
throughout the orbit due to the inclination angle, or the direction of motion, and most certainly they could not have 
known in which quadrant the ascending node would be. It is quite likely, however, that since the period was known the 
measurements are not entirely independent, and some may well be wild guesses by inexperienced observers. We note, 
for example, that the quadrants of the individual observations as published are all consistent with the modern orbit, 
with the exception of two that need to be reversed. Rather than implying that the observers were able to determine 
the correct orientation in every case by detecting the slight difference in brightness, it seems more likely that they 
adjusted their 6 values by ±180° prior to publication, as needed for consistency with the initial measurements. This 
practice is still followed today, but indicates at the very least that the Greenwich observers were paying close attention 
to the agreement with the spectroscopic orbit. Three observers contributed almost 85% of all measures shown in 
Figure [TTI About half of them were carried out by W. W. Bryant, and another third by W. Bowyer and T. Lewis. No 
systematic differences are seen among them. The very detailed observing notes they left indicate that they also made 
similar measurements of other very close binaries on some of the same nights on which they observed Capella, and to 
our knowledge those results have not been called into question. 

After more than a century the Greenwich observations remain a mystery. H. F. Newall, who in his earlier spectro- 
scopic study at that Observa tory had disc overed the binary nature of Capella and encouraged the visual observers 
to turn their attention to it (jNewalll [19000 . later expressed some doubts a bout those m easurements in meetings of 
the Royal Astronomical Society, but believed that at least some were real (iNewalll llQO^. Bryant himself, the most 
active observer, took a similar view. In his own words (Bryant 1907), "though many of the observations published are 
probably worthless, and may be ascribed to optical or psychological causes and to want of experience in judging of 
the suitability of observing conditions for what was practically a unique observation, yet there is little doubt that a 
certain proportion of the observations can be regarded as genuine." 

NOTES ON THE ASTROMETRIC OBSERVATIONS 

We describe here the particulars of the astrometric measurements for Capella, including corrections to some of the 
original sources, or re asons t hat led us to reject some observations from th e final solution. 

\Andersor\ Hi 9 2(1 ) and \MerTil i (1922) : The six measurements reported by lAndersonI ()1920f) were re-reduced and are 
superseded in the paper bv iMerrilll (1922). The P. A. for first of these six (on JD 2,422,323.65) is admittedly very 
rough, and has been assigned an uncertainty of 10°, following the recommendation of the authors. The last of those 6 
measurements (on JD 2,422,438.63) was made during full daylight, and thus an angular separation was not reported 
because it was affected by the blue backgro und tha t shift ed the effective wavelength. The P.A. should be unaffected. 
The measurements reported in Table II by IMerrilll (|1922t l are nightly averages, but smearing should be insignificant 
given that the observing intervals were always less than 3 h ours. The relative uncertainties used here for all these 
measurements account for the weights given W iMerrilH ()1922f ). and are scaled as described in §[4l We note in passing 
that the dates o f the se observation s have been translated into Bessellian epochs incorrectly in the orbital studies of 
iMcAlisteil (|1981| ) and lBarlow et al.l ()1993l ). The values used by those authors are exactly 0.5 days too small, possibly 
as a result of overlooking the fact that prior to 1925 astronomical time (GMT) was counted from noon, rather than 
midnight. The change in the P.A. and separation of Capella over half a day can be up to 2?4 and 0.45 mas, respectively. 
\KulaaiT\ IT97G) : The original posit ion angle for th e interferometric observation on 1968 February 18/19, listed as 
303?5± 1?0, is apparently incorrect ()McAlisteij|1981h . and was changed by the latter author to 313?5± 1?0 following a 
private communication from W. S. Finsen. Even after this change, the P.A. measurement gives a very large residual, 
and is exclud e d from our solution. Nothing wrong is seen with the separation measurement on this date. 
\Gezari et al\ 1197&): A speckle measurement of the angular sepa ration is g iven with no epoch. It is compared in 
Table 1 of the original paper with predictions from the orbit by M errilll ()1922l ) for JD 2,441,034.0, which may then be 
taken as the date of observation. This date is likely to be rounded off, however, since the star was below the Palomar 
Observatory local horizon at the time. It corresponds to 1971 March 23, which also disagrees with indications in the 
text saying that "obse rvations were niade d uring twilight on 11 dates in 1971 April, June, and October". A later paper 
by the same authors ( Labevrie et al.|[l974 see below) does report measurements from the same site on 1971 March 
23, so it seems likely they correspond to the same observations, although the angular separation reported is slightly 
differe nt (possibly as a resu lt of new reductions). We assume here that the latter measurements supersede the original 
one bvlGezari et all (Il97l) . 

\Labevrie et al\ U97A ): The dates for the ten speckle observations given here are only reported to 0.01 of a y ear, which 
is insufficient given that orbital motion is ~10° and ~2.1 mas over this time. Previous investigators (e.g.. iMcAlisted 
[l981; Barlow et al. 1993) have used th e dates at face value, n ot surprisingly with poor results. Furthermore, the first 
seven measurements were reported bv iLabevrie et ahl ()1974[ ) with precisely the same date (1971.23, presumably in 
Bessellian years), and were therefore simply averaged by subsequent investigators. It is clear from the text that the 
first four of these were taken on the night of 1971 March 23, probably during twilight at Palomar since the star was 
setting. We assign to these observations the date JD 2,441,034.67. This is unlikely to be in error by more than 2 hours, 
which is good enough for our purposes. Inspection of measurements for other binaries from this source, as well as those 

Both are by the same observer (T. Lewis). Interestingly, on one date the quadrant was stated in the observing log as being ambiguous, 
and on the other the measurement was qualified as uncertain. 
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of Capella, indicates that the position angles are systematically too small by several tens of degrees, a fact also noticed 
by others (despi te the date problem). This may be due to a calibration error in these early speck le observations 
( Mc Alist eil 1 1 9 8 If ) ■ and we have chosen not to use those angles here. The remaining observations of iLabevrie et al] 



19741) were l ikely ta ken a few nights later, but sufficiently precise dates are unavailable so we have excluded them. 



\Blazit et al\ \ 197 7 A) : The results of these interferometric measurements are taken from Figure 1. The date for the 
last one (1977.207) is given to one more significant digit in the text (1977.2066). Dates are assumed to be Bessellian 
years, more common in binary star publications, and also because the assumption of Julian years for the 1977.2066 
measurement places the star below the horizon. All dates are converted to Julian days for use in the present work. 
The observation on 1977.125, which corresponds formally to JD 2,443,189.6, is considered to have been made 0.1 days 
earlier, when the star was at a more plausible hour angle. This is well within the uncertainty in the original date 
(0.001 years), which corresponds to ±0.36 days. We note that the change in 6 over 0.36 days is about 1?6 on this 
date (p only changes by 0.2 mas), so our date correction is not entirely insignificant. Uncertainties for the 9 and p 
measureni e nts ar e estimated graphically from Figure 2 of this source. 

\McAliste^ il97ll ): The position angles of these speckle measurements were later slightly adjusted by the original 
author. 

I Weiaeli IW/A ): Capella was unresolved in this speckle measurement made at an unspecified date with the 1.8m 
telescope at the Asiago Obse rvatory (Italy). 

\McAlister & HendrA i ' 198 2d ): The first of these speckle measurements, which were made at the resolution limit of the 



2.1m telescope at Kitt Peak Observatory, gives large residuals both in p and 9, and was excluded from the solution. 
\Hege et all ^198 A ): The original speckle observations gi ving p = O " 0A2 ± y 001 and 9 ^ 331° ± 2° (quadrant reversed 
for the present work) were reduced independently by iBagnuolol ()1983af ). who obtained p — 0'.'0420 ± O'.'OOIS and 
9 = 333?0± 1?2. We have adopted here the average of these determinations, with the more conservative uncertainties. 
The date assigned to this observation (not precisely given by either author) is JD 2,444,637.6, which corresponds to a 
reasonable hou r angle at the Kitt Peak Observatory during twilight on UT 1981 February 2. 

\Koechlin et al! {198a): The first three interferometric measurements are apparently repeated from an earlier report 
by Kocchlin et al. (1979), although they are slightly different in two cases. We consider them to supersede the original 
measurements. The date for the last of these three is also different from that given in the original, and is the one used 
here. 

il98't) : These speckle measurements were later slightly corrected by the original authors. 
198^ ): The speckle observation on JD 2,446,895.62 (Bessehian epoch 1987.2717) gives very large 
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residuals in both and p, and is excluded from our solution. 
\Isobe et~ai\ H99(\ ): Capella was unresolved in this speckle measurement made on 1988.8031 (Bessellian epoch) at the 
2.12m telescop e at th e San Pedro Martir Observatory (Mexico). 

\Hummel et al\ 1(1994 ): Dates given in Table 2 of this source have been converted from Bessellian years to Julian dates 
for use in the present work. The interferometric visibilities on which the results are based have not been published, 
but wer e condensed into and published as 9 and p measurements for each night of observation (at 8*^ UT). 
\Baldwin et all Iil99d) : The P. A. and separation values used in the present work were determined graphically from 
Figure 2 of this source, since they were not given explicitly. The epoch of observation for the first interferometric 
measurement is deduced from the calendar date mentioned in the text, along with the time during the night (2^ UT) 
from Figure 3. The second measurement is assumed to be taken at a similar time during the night of 1995 September 
29^ 

\Youno & Duvre^ 1^200 A ): Exact times of observation for these direct-imaging HST measurements were not given in 
the original paper, but have been recovered from the HST online archives. The position angles are off by nearly 40° 
for unknown reaso ns, and have not been used here. 

\Kraus et al\ \200A ): The time of observation is taken to be the middle of the interval given for each measurement in 
Table 2 of this work. 

CORONAL ABUNDANCES IN CAPELLA 

There is a significant body of literature on X-ray observations of Capella, which is not surprising given that it 
is the brightest steady coronal source in the X-ray sky. The available abundance determinations made from coro- 
nal emission-li ne fluxes have been collect ed in Table [20l and have been adjusted to conform to the standard solar 
abundances bv iGrevesse fc Sauvall (|1998f l. for consistency with the photospheric determinations (see §[7|). Coronal 
abundances in the Sun are known to depend on the first ionization potential (FIP), although the large-scale element 
fractionation processes presumably involved are not completely understood. Elements having low FIP (< 10 eV) are 
overabundant compared to those with high FIP (> 10 eV). We indicate the FIP values in Table \20\ along with the 
average abundance for each element in the last column. The FIP trend has also been seen in some active stars, while 
others appear to show the opposite effect. F or Ca pella there ha s been some debate about this issue , som e authors 
claiming a mild FIP effect d Brickhouse e t al.ll2001l : rAudard et ai]l2003t lArgiroffie^aH 120031: Iciiet al.ll2006j V. or even 



an inverse effect |Bauer fc B rcgman 1996), and others claiming no effect at all (jFavata et al.Hl997tlAudard et al.ll2001l : 



iMewe et al.ll2001[ ). Some of this disagreement may stem from the use of different solar abundances, which is why we 
have homogenized all measurements here. As indicated in §[7l most authors agree that nitrogen is overabundant. The 
average abundance of the remaining high-FIP elements (C, O, Ar, Ne) is —0.30 ± 0.04, while that of the low-FIP 

Some systematic differences between studies are apparent, but tlieir investigation is beyond the scope of the present work. 



Capella 



27 



0.2 - 



X 



0.0 - 



0.2 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Ca T 



Ni,Mg, 
Fe,Si 



N 



S 



0.4 - 



C 



.0 
i 4 



Ar 



Ne 



Photospheric [m/H] 



_] I I L. 



10 15 20 

First ionization potential (eV) 



25 



Fig. 18. — Measurements of the coronal abundances in Capella from X-ray observations (combined flux of the two components) as a 
function of the first ionization potential of each element. The values shown are averages of several sources given in Table [20l and are taken 
from the last column of that table. Error bars represent the error of the mean, and all measurements are referred to the solar abundances 
as listed by Grevessc & Sauval (1998). Nitrogen is marked with a different symbol because it is expected to be enhanced due to convective 
mixing in the pri mary star occurring at the first dredge-up. The overall photospheric composition of Capella, [m/H] = —0.34 ± 0.07 
l)McWilliamlll99"0t ). is indicated with the dashed line and gray error box, and refers to the same solar abundances. 

elements (Ca, Ni, Mg, Fe, Si) is -f0.01±0.03, although it is likely that the true uncertainty is larger than indicated by 
these statistical errors. Sulfur has been excluded since its first ionization potential is close to the conventional dividing 
line of 10 eV. The FIP effect is thus clear in Capella: low-FlP elements are enhanced by a factor of ~2 relative to 
high-FIP elements (see Figure ITSl) . Furthermore, the ev idence in the Sun indicates that high-F IP elements show better 
average agreement with the photospheric composition ()Meveiill993l : iFeldman fc Widindl2007l . and references therein). 
If the same were true in Cape lla, then the coronal abundance as measured from the high-FIP elements supports the 
photospheric determination of iMcWilliamI (jT990l ). which gives an overall metallicity of [m/H] = — 0.34± 0.07. Under 
the same assumption, the enhancement of N relative to the high-FIP (or photospheric) composition would also be 
about a factor of 2, or slightly higher. Considering that this measure is diluted to some degree by the flux contribution 
fro m th e (presumably normal composition) secondary, it appears to agree at least qualitatively with model predictions 
bv iBer telli et al. (2008), which indicate a N enhancement of a factor of ~3.1 for a post first dredge-up star. The 
average C/N ratio from Table [201 (note that since both C and N have high FIP, any residual FIP effect is more likely 
to cancel out) is approximately 0.52, though with a large uncerta inty of at least . 2-0.3 due to the scatter in C. An 
independent and perh aps more accurate es timate was reported bv lSchmitt fc Nes"i ()2002D as C/N = 0.91 ± 0.04. The 
ratios predicted by the I Bertelli et al.l (|2008f) models are 0.67 after first dredge-up and 3.27 for the unevolved secondary, 
again roughly consistent with the measured values for the combined flux. Other models considered in ii l9.1l give similar 
values. 

Although a somewhat clearer and consistent picture than one might glean from a casual look at the literature 
seems to emerge from the above regarding coronal abundances in Capella, uncertainties are still rather large and our 
understanding of heating and element separation in the outer layers of stellar atmospheres is still limited. 
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TABLE 1 

CfA radial velocity measurements of Capella. 



HJD 




RVa 


RVb 


CTA 




(O-C)a 


(O-C)b 




(2,400,000+) 


Julian Year 


(kms"-'-) 


(kms"-"-) 


(kms^-"-) 


(kms^-"-) 


(kms^-') 


(kms-i) 


Phase 


46475.6179 .... 


1986.121 


47.98 


8.93 


0.61 


1.21 


-0.42 


-1.53 


0.8781 


46494.5383 .... 


1986.173 


53.18 


5.27 


0.68 


1.37 


-0.65 


+0.30 


0.0600 


50388.7228 .... 


1996.834 


3.65 


55.02 


0.54 


1.08 


-0.01 


-0.62 


0.4963 


50391.7233 .... 


1996.843 


4.42 


56.57 


0.43 


0.86 


+0.45 


+1.25 


0.5251 


50438.8498 .... 


1996.972 


55.98 


4.68 


0.56 


1.12 


+0.57 


+ 1.31 


0.9781 



Note. — Tabic [T] is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



TABLE 2 

Spectroscopic orbital solutions from published velocities for Capella. 





7 


Prim/Sec offset 


Ka 


Kb 


T cr\ (TB 


Time span 


Source 


(kms-i) 


(kms-l) 


(kms-l) 


(kms-l) 


{HJD-2,400,000) (kms-i) (kms"!) A^a Nb 


(days) 



+30.15 


± 


0.14 




25.89 


± 0.18 








14865.54 


± 


0.14 


0.77 




31 




1486.0 


+27.19 


± 


0.55 




26.89 


+ 0.83 








1 A nan q a 
14969. o4 


± 


0.48 


2.20 




23 




121.0 


-hoU.D / 




n 07 
U.z / 


U 


zo.oo 






_|_ 
it 


1 /I 


i / 4DD.0U 


± 


U.zo 


l.oD 


o.4i 


oO 




1 nn/i Q 
1UU4.0 


+29.66 


± 


0.08 





27.06 


± 0.20 


34.34 


± 


1.49 


23085.33 


± 


0.14 


0.27 


3.36 


7 


4 


55.9 


+31.44 


± 


0.17 





26.02 


± 0.24 


31.09 


± 


0.82 


28180.59 


± 


0.15 


1.01 


3.10 


36 


16 


753.7 


+29.30 


± 


0.12 




26.46 


± 0.17 








31821.46 


± 


0.12 


0.83 




45 




8142.9 












53.23 


± 0.63 






34214.36 


± 


0.50 




2.42 




19 


802.9 


+29.50 


± 


0.15 




26.10 


± 0.20 








42119.27 


± 


0.15 


0.64 




18 




274.0 


+29.57 


± 


0.14 


+0.84 ± 0.45 


26.35 


± 0.19 


29.75 


± 


0.59 


44927.73 


± 


0.11 


0.86 


2.58 


42 


39 


1534.8 


+29.87 


± 


0.13 





26.14 


± 0.18 


30.43 


± 


0.50 


45031.92 


± 


0.12 


1.33 


3.85 


105 


91 


1902.9 


+31.66 


± 


0.59 




25.48 


± 0.83 








46593.65 


± 


0.51 


2.92 




29 




957.6 


+28.79 


± 


0.87 




23.0 


± 1.1 








46594.46 


± 


0.86 


3.62 




24 




957.6 


+30.50 


± 


0.12 





26.29 


± 0.17 


27.06 


± 


0.46 


42119.37 


± 


0.11 


0.46 


1.31 


15 


15 


211.7 


+29.14 


± 


0.10 





25.96 


± 0.11 


27.70 


± 


0.41 


44615.96 


± 


0.12 


0.47 


1.38 


23 


14 


4532.6 


+30.36 


± 


0.47 


-6.44 ± 0.95 


22.43 


± 0.61 


35.83 


± 


1.21 


48568.57 


± 


0.34 


2.77 


2.64 


40 


16 


1562.5 


+29.66 


± 


0.04 


+0.27 ± 0.08 


25.96 


± 0.04 


26.27 


± 


0.09 


50857.21 


± 


0.03 


0.46 


0.91 


162 


162 


5037.2 



Campbell (19011 

Newall (190qf~^ 

Goos (1908) 

^nford (1922) 

Struvc (1939) 

Struvc & KUbv (1953 )" 

Wright (19^^ 

Batten & Erceg (1975) '^ 

^ficn et al. (1985)'* 

Beavers & Bitter (1986) 

'Sficficrbak ov et al. 0990) ° . . . 
~hclierbakov et al. (1990) ^ . . . 

Batten et al. (1991) 

Barlow et al. (1993) 

j<atsova_^_Sclierbakov Q998) 
CfA observations (this work) . 



Note. — The follow ing radial-velocity measurements were rejected in com puting these solutions: JD 2,445,002.6506 jShen e^^y ^85 |. primary): 2 ,444.930.825. 2,445,366.634 
JBeavers fc Eitte jIToSel . primary); 2,445,315.781, 2,445,315.793, 2,445,320. 718 j Bcavcrs fc Eittcr" 1986|, secondary); 2.446.345.71 llShcherbakov et al."1990'. primary from Ho I A10830 
measurements); 2,449,242.545, 2,449,315.281, 2,449,316.458, 2,449,316.480 l IKats ova & Schcrbakov 199j, primary); and 2.449.443.384 lIKatsova fc S chcrbakov 1998, secondary). 
" The date of first of the observations on 1951 January 2 has been corrected, following (Heint^ fll975D The measurements were published only in the form of velocity differences, and 



dates were given modulo the orbital period (original dates not available). ° The velocities measured by these authors were superseded by those of lBatten et all (I1991D ■ who re-measured 
most of the original photographic plates taken bv lBatten &: Erceei jl975D .^ We consider the velocities measured by these authors to be superseded by those reported bv lBeavers fc Eitterl 



who added one more year of observations.^ The elements reported here are based on velocities measured from photospheric features in Capella. 
are based on velocities measured from the He 1 A10830 line, which is not photospheric. 



The elements reported here 
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►a 

CD 
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TABLE 3 

Historical radial velocity measurements of the primary component of 

Capella. 



JD 




RVa 


O-A 


(O-C)a 






(2,400,000+) 


Julian Year 


(km s^-"-) 


(km s"-"-) 


(km s^-"-) 


Phase 


Source 


13804.0360 .... 


1896.671 


36.50 


0.75 


-0.72 


0.7939 


1 


13820.0050 .... 


1896.715 


53.90 


0.75 


-0.84 


0.9474 


1 


13837.0030 .... 


1896.761 


50.40 


0.75 


+0.31 


0.1108 


1 


13839.0290 .... 


1896.767 


47.00 


0.75 


-0.90 


0.1303 


1 


13876.8650 .... 


1896.870 


4.30 


0.75 


+0.15 


0.4940 


1 



Referenc es. — (1) [C ampbell l ilQOlD ; (2) Ncwall (190p) ; (3) IGooj (I1908D: (4) IStruve 
{1939); (5) , Struve fc Kilby 11953i): f 6) [Beave rs & EittciT lligSeD : (7) IShcherbakov et al, 
n99(]|) : (Sl lBatten et al.l | |199]1) : l'9) IBarlow et al.l ||1993I'I . 

Note. — Table[3]is presented in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 



TABLE 4 

Interferometric measurements of Capella published in polar coordinates. 



JD 9 p ag ap (O-C)g (O-C)p 

(2,400,000+) Julian Year (°) (mas) (°) (mas) (°) (mas) Phase Source 



22323.65 1919.996 148.0 42.8 10.00 0.54 -8.66 +1.16 0.696 1 

22368.63 1920.119 5.0 46.9 1.73 0.80 +1.46 -0.35 0.129 1 

22369.65 1920.122 1.0 46.2 1.44 0.67 +0.94 -0.16 0.138 1 

22370.66 1920.125 356.4 45.4 1.44 0.67 -0.08 -0.09 0.148 1 

22399.63 1920.204 242.0 51.7 1.73 0.80 +0.91 +0.00 0.427 1 



References. — dl lMerrilil U922f) : (2 HKulaginl 1119701) : (•3') IBlazit et al.l J1977all : (41 IHummel et al.l U994H 
Note. — Table^is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion 



is shown here for guidance regarding its form and content. 



TABLE 5 

Interferometric measurements of Capella published in Cartesian coordinates. 



JD 




Ax 


Ay 






{0-C)Aa: 


iO-C)Ay 






(2,400,000+) 


Julian Year 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


Phase 


Source 


43773.6980 .... 


1978.723 


+40.0 


+19.4 


6.70 


0.76 


-8.07 


-0.49 


0.904 


1 


43792.6460 .... 


1978.775 


+14.3 


+50.8 


6.70 


1.26 


-1.09 


+0.21 


0.086 


1 


43798.6290 .... 


1978.792 


-2.8 


+47.9 


4.02 


1.26 


-1.00 


+0.17 


0.143 


1 


43805.6100 .... 


1978.811 


-23.8 


+36.4 


2.68 


0.51 


-2.39 


-0.28 


0.210 


1 


43813.5880 .... 


1978.832 


-34.2 


+16.8 


5.36 


1.26 


+4.85 


+0.48 


0.287 


1 



Ref erences. — (1) iKoechlin eTall ||1979D : (2) iKoechlin eTall ||T983); (3) iBaldwin erall [1199611 : (4) iKraus et all 
(f2005h . 

Note. — Table [5] is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



Capella 



TABLE 6 

Speckle and imaging measurements of Capella. 



JD 




e 


P 


0-9 




{0-C)g 


(o-c)p 






(2,400,000+) 


Julian Year 


n 


(mas) 


n 


(mas) 


n 


(mas) 


Phase 


Source 


41034.6700 .... 


1971.224 




53.00 




1.36 




-0.85 


0.572 


1 


41034.6700 .... 


1971.224 




54.00 




1.36 




+0.15 


0.572 


1 


41034.6700 .... 


1971.224 




55.00 




1.36 




+1.15 


0.572 


1 


41034.6700 .... 


1971.224 




56.00 




1.36 




+2.15 


0.572 


1 


42643.1500 .... 


1975.628 


35.00 


53.00 


5.76 


3.40 


+3.75 


-2.78 


0.035 


1 



Ref erences. — (1) ILabeyrie et al.| Jl974|): (2) IB lazit et aTl lll977bD; (3) |McAlistei| ||1977I): f4 ) IMcAlisteil 
a978l): (5) [M cAlistcr fc DcGioia (1979); (6) Balcga lTTiknonovl II1977D: (71 McAlistcr & Fckcl |'198 q); (8) 
IMorgan etEO^J ^gSOV (9) McAlistcr & Hcndrv (1982a'): f 10') IMcAlister fc Hend ry (1982b): (11) McAlistcr et al.l 
ai983l): (12TlDudinov et al.l (11982): (13) Hcgc ct al. (1983) and'BaEnuolrf ('1983a): (14) BalcEa et al. (1984); (15) 
IMcAlister et al.l 11198 ' H); (16 )|BalcEa & Balcga (1985); (17) Hartkopf ct al. (2000): (18) Balcga & Balcga (1987); 
(19) 'Bonncau ct al.' (19861; (20) McAlistcr ct al. (1989); (21) McAlistcr et al.l (119931) ; (22) iMcAlistor ct aO 
11990); (23) Hartkopf ct al. (1992); (24) Horch ct al. (1999); (25) Y oung &: Dupred ( 12003) (direct imaging rather 
than speckle); f26) ,Scardia et al.l (1200Q) ; (27) iBalcEa ct al.i (i2004i) . 

Note. — Tablc|6]is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion 
is shown here for guidance regarding its form and content. 



TABLE 7 

Hipparcos MEASUREMENTS OF CAPELLA AND CORRESPONDING O — C 
RESIDUALS. 



HJD 




Av 


0"Au 


O-C 




(2,400,000+) 


Julian Year 


(mas) 


(mas) 


(mas) 


Phase 


47921.1897 .... 


1990.079 


-6.53 


3.54 


-7.11 


0.7749 


47921.2922 .... 


1990.079 


-0.61 


2.57 


-1.19 


0.7759 


47933.1656 .... 


1990.111 


-2.90 


1.85 


-0.99 


0.8901 


47933.2030 .... 


1990.111 


-1.76 


3.50 


+0.16 


0.8904 


47965.5793 .... 


1990.200 


-0.71 


1.68 


+0.87 


0.2017 



Note. — Tabled is presented in its entirety in the electronic edition 
of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



TABLE 8 

Summary of data sets used in the orbital solution for Capella, along with 

TYPICAL uncertainties. 



Source 



Measurements 



Spectroscopic observations 

Campbell ( 1901) 

Newall (lOOoT"*^ 

Goos (1908)7 

^ruve (1^939) 

Struve k K ilbv (1953) . . . 

Beavers fc E itter ( 1986) . . 

aicherbakov ct ah (1990) 

fatten ct al. ( 19M) 

^arl^^_et_al^ 

^^^^thislpaper] 
Astrometric observations 

Merrill (1922) 

KulaginTigrO) 

Blazit et ar7l977a) 

Koechlin et al. (1979, 1983) 

Hummel ct al. ( 1994) 

Baldwin et al. (1996) 

Hipparcos (Pcrrvman et al. 

Ypuns fc Duprce (2002 ) 

Kraus ct al. (200^ 

Speckle observers 



1997) 



31 measures of -RVa 







74 




23 measures of -RVa 




2 


12 




35 measures of -RVa 




1 


52 




36 measures of -RVa 







97 




45 measures of -RVa 







91 




105 measures of -RVa 




1 


27 




29 measures of -RVa 




3 


14 




15 measures of -RVa 







48 




23 measures of -RVa 







45 




162 measures of -RVa and -RVb 


0.44 


, 


89 


16 measures of 8 and 15 of 


P 


1.50 


, 


62 


11 measures of 8 and 12 of 


P 


1.10 


, 1 


52 


3 measures of 8 and p 




5.98 


, 3 


39 


25 measures of Ax and Ay 




2.67 


, 


67 


40 measures of 8 and p 




0.13 


, 


11 


2 measures of Ax and Ay 




0.93 


, 1 


57 


41 measures of Ad 




2 


82 




4 measures of p 




2 


19 




5 measures of Ax and Ay 




0.26 


, 


19 


47 measures of 8 and 53 of 


P 


1.18 


, 1 


19 



^ RMS residual of unit weight from our global fit, in the natural units of each type of observation 
{kms^^ for RVa and RVb, degrees for 9, and mas for p, Ax, Ay, and Av). 
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TABLE 9 
Orbital solution for Capella. 



Parameter Value 



Adjusted quantities 

P (days) 104.02173 ± 0.00022 

a" (mas) 56.408 ± 0.025 

i (deg) 137.212 ± 0.051 

f^.72000 (deg) 40.421 ± 0.064 

T (HJD-2,400,000) 47528.514 ± 0.016 

7 (kms-i) +29.653 ± 0.035 

Ka (kms-i) 26.005 ± 0.036 

Xb (kms-i) 26.260 ± 0.087 

''phot i^^) 2.08 ± 0.69 

Aa* (mas) -0.26 ± 0.80 

AS (mas) -0.30 ± 0.57 

A^i* (mas yr-i) +0.32 ± 0.99 

Ans (mas yr"!) -0.08 ± 0.59 

fp'^ 1.0400 ± 0.0035 

Aab (kms-i)'= +0.267 ± 0.079 

Ai (CfA - Campbell ^901) ^ (kms'M -0.48 + 0.14 

A2 (CfA - Ncwall (1900)) (kms-^) +2.12 ± 0.45 

A3 (CfA - Coos (1908) ) (kms-i) -0.93 ± 0.26 

A4 (CfA - Struve (1939)) (kms"!) -1.82 ± 0.17 

A5 (CfA - Stru ~ fc Kilbv (1953 )) (kms-^) +0.34 ± 0.14 

Ae (CfA - Beavers fc Eitter (1986)) (kms-^) -0.37 ± 0.13 

A7 (CfA - Shch erbakov et al. (1990)} (kms-^) . . . -2.37 ± 0.59 

Ag (CfA - Batt en et al. (1991) ) (kms"!) -0.90 ± 0.13 

A9 (CfA - Barlow et al. (1993) ) (kms"!) +0.50 ± 0.10 

Derived quantities 

Ma (Mq) 2.466 ± 0.018 

Mb (Mq) 2.443 ± 0.013 

q = Mb/Ma 0.9903 ± 0.0036 

a (10'5 km) 110.06 ± 0.23 

a (AU) 0.7357 ± 0.0015 

TTorb (mas) 76.67 ± 0.17 

Distance (pc) 13.042 ± 0.028 

/X* (mas yr-i) +75.84 ± 0.99 

lis (mas yr-i) -427.21 ± 0.59 

(^B/^A)ffp 1.171 ± 0.057 



^ Scale factor for the angular separation measurements by IMerrilll jl922D and iKulaginI 
jl970D .^ Primary/secondary offset for the CfA velocities, in the sense (primary minus 
secondary) . 
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TABLE 10 

Relative brightness measurements for Capella. 



Source 



Youn g & Duprcc (2002) 

Youn g & Duprcc (2002) 

Bagnuolo fc S owcU (1988)'' 

Griffin fc Griffi n (1986) 

Hummel ct al. ( 1994) 

Bagnuolo & SowcU Tl988) =" 

CfA spectroscopy (this paper) 

Hege et a l. (1983)^^ 

Hipparcos (this paper) 

_Bagnuolo & SowcU (1988) =" 

Bagnuolo & SowcU (WW)'^ 

Bagnuolo & SowcU (1988)^ 

_B^nuolo & Sowcll (1988)=" 

Hum mel ct al. ( 1994) 

Frost & Rushforth (1979, priv. comm.)'* . . . 

Bagnuolo (1983a) 

Bagnuolo (1983b) 

Griffin fc Griffin (1986) 

Strassmeier fc Fekel (1990) 

Hummel et al. (1994) 

Bald win et al. (1996) 

You ng (1999)" 

Kra us ct al. ( 2005) 

^^^cncdctto^_^onneau (199^) 



Wavelength Bandwidth 



(^m) 


(pm) 


■^cool /'' 


?hot 


0.1400 


0.0200 


0.059 


± 


0.012 


0.2750 


0.0500 


0.25 


± 


0.04 


0.4108 


o nooQ 


0.61 


± 


0.02 


0.4425 


0.0988 


0.69 


± 


0.05'' 


0.4500 


0.0200 


0.77 


± 


0.07 




n 0909 


0.82 


± 


0.02 


U.Oloo 




0.677 


± 


0.023 


^200 


0.0100 


0.64 


± 


0.03 


0.5275 


0.2378 


0.854 


± 


0.041 


0.5478 


0.0253 


0.93 


± 


0.02 


0.5478 


0.0253 


0.87 


± 


0.02 


0.5478 


0.0253 


0.92 


± 


0.02 


0.5478 


0.0253 


0.91 


± 


0.02 


0.5483 


0.0200 


0.87 


± 


0.04 


0.55? 




0.65 


0.5500 


0.0100 


0.82 ± 


0.02 


0.5520 


0.0200 


0.89 


± 


0.01 


0.5544 


0.0898 


0.87 ± 


0.05^= 


0.6430 


0.0084 


0.91 ± 


0.05^= 


0.8103 


0.0200 


1.05 ± 


0.05 


0.8300 


0.0400 


1.11 ± 


0.05*= 


1.3 


0.16 


1.38 ± 


0.06 


1.65 


0.3 


1.49 ± 


0.10 


2.15 


0.11 


1.54 ± 


0.35 



=" Measurements have not been transformed to the standard Stromgren system.'' Un- 
certainty not originally reported, and assigne d arbitrarily for th is work.'^ Supersedes 
the result of ICocke et al.l (|1983l ).'' Cited bv iHeee et"al] (| 19831 ). The passband is 
uncertain. 



TABLE 11 

Limb-darkened angular diameter measurements for Capella. 



Source^" 


eA 

■^LD 
(mas) 


"-^LD 
(mas) 


Di Benedetto & Bonneau (1991) . . . 

Hummel et al. (1994) 

Kraus et al. (2005) 


9.5 ± 2.3 
8.46 ± 0.10 
9.14 ± 0.62 


6.18 ± 0.42 
6.32 ± 0.30 
5.92 ± 0.82 


Weighted average 


8.47 ± 0.40 


6.24 ± 0.23 



^ The original measurements have been homogen ized by uni- 
formly applying limb-darkening coefBcients by Ivan Hammd 
(|l99a i. 
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TABLE 12 

Summary of photospheric abundance determinations for 
Capella. 



Element or ratio 


- 

Primary 


Secondary 


Source 


[Si/H] 


-0.32 (3) 




1 


[Ca/H] 


-0.26 (2) 




1 


fSc/Hl 


-0.73 (1) 




1 


[Ti/H] 


-0.28 (8) 




1 


[V/H] 


-0.09 (2) 




1 


Fc /HI 


—0.20 ± 0.22 




1 


[Co/H] 


-0.02 (2) 




1 


[Ni/H] 


-0.63 (7) 




1 


[Sr/H] 


+1.18: (1) 




1 


[Y/H] 


-0.20 (1) 




1 


[Zr/H] 


-0.77 (1) 




1 


[La/H] 


-0.32 (1) 




1 


[Nd/H] 


-0.27 (1) 




1 


[m/H]*^ 


-0.34 ±0.07 




1 


logAr(Li) 


1.3 ±0.2 


3.2 ±0.3 


2 


12C/13C 


27 ±4 




3 


C/N'^ 


0.57 ± 0.06 


3.30 ±0.16 


4 



Refe rences. — (l)IMcWilliaml lUgga). converted to the abundance 
scale of lGreves sc fc SauvalMl99^: f 2') |Pilachows ki fc SowcT? fT99^ , re- 



computed here to account for temp erat ure and g ravity differences, and 
NLTE effects. Q^ iTomkin et alj U97^ : (4) [Linskv ct al.i a995i) : 
^ Values in parentheses represent the number of spectral lines measured 
bv iMcWilliamI jl990D . AH abundances relative to hydrogen have been 
converted to the scale of lGrevesse fc Sauvaj < |199^ Average metal- 
licity excluding [Sr/H].^ Ratios derived from transition layer lines, be- 
lieved to be representative of the photospheric values (see text). Un- 
certainties for the line fluxes of CIVA1550.8 and N V A1238.8 us ed to 
compute the C/N ratios were not reported bv lLinskv et all (119951) sep- 
arately for the two components; we have somewhat arbitrarily assumed 
here the same relative uncertainty as the ratio calculated from the total 
flux of each line i jLinskv et al.iri995l . Table 3), and then doubled that 
of the primary on account of its fainter lines (sec §[TJ- 
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TABLE 13 

Apparent magnitudes, color indices, and inferred effective temperatures for Capella. 



Passband 



^cool / '-hot 



Combined light'' 



Primary star 



Secondary star 



Apparent magnitudes 




















a (magj 


n 71 ^ ^ n C\^^ ( A\ 


u.o / u 


-1- 


U.UlO ) 


1 S9n 


-1- 


u.uoo 


1 A^fi -t- 


ft 041 


V t^mag; 


n + n ni 1 (q\ 

u.ooo zn. u.u-Li- } 


0.073 


-|- 




0.892 


-|- 


0.016 


0.763 it 


0.015 


He (mag) 


0.91 ± 0.05 (1) 


-0.346 


± 


0.031 (3) 


0.459 


± 


0.044 


0.357 ± 


0.042 


Ic (mag) 


1.080 ± 0.035 (2) 


-0.734 


± 


0.018 (3) 


-0.022 


± 


0.025 


0.061 ± 


0.026 


■h (mag) 


1.38 ± 0.06 (1) 


-1.337 


1 


U.UyU [2) 


—0.745 


1 


0.092 


—0.396 it 


0.094 


H2 (mag) 


1.49 ± 0.10 (1) 


-1.57 




0.18 (1) 


— 1.02 




0.19 


—0.58 it 


0.19 


K2 (mag) 


1.54 ± 0.35 (1) 


-1.795 


± 


0.082 (2) 


— 1.25 


± 


0.14 


—0.78 it 


0.17 


Color indices 




















B-V (mag) 




0.797 




U.UZO 






u.uo / 


0.693 it 


0.044 


V-Rc (mag) 




0.419 


-1- 








n C\A7 


0.406 it 


0.045 


V-Ic (mag) 




0.807 


-1- 


u.uzo 






u.uou 


0.702 it 


0.030 


(R-I)c (mag) 




0.388 


± 


0.036 


0.481 


± 


0.051 


0.296 it 


0.049 


V-J2 (mag) 




1.410 


± 


0.091 


1.637 


± 


0.093 


1.159 it 


0.095 


V-H2 (mag) 




1.65 


1 

± 


0.18 


1.91 


1 

± 


0.19 


1.35 it 


0.19 


V -K2 (mag) 




1.868 


± 


0.083 


2.14 


± 


0.14 


1.55 it 


0.18 


Temperatures (K) inferred from the color indices, assuming [m/H 


= 


-0.34 ±0.07 












B-V 




5138 


± 


76 


4862 


± 


124 


5391 it 


128 


V-Rc 




5219 


± 


199 


5142 


± 


258 


5294 it 


263 


V-Ic 




5224 


± 


81 


4927 


± 


85 


5582 it 


120 


(.R-I)c 




5238 


± 


94== 


4748 


± 


195 


out of range 


V-J2 




5244 


± 


159 


4861 


± 


144 


5806 it 


119<= 


V-H2 




5423 


± 


254 


5092 


± 


223 


5890 it 


232 


V-K2 




5212 


± 


105 


4905 


± 


146 


5657 it 


244 


Weighted average 




5207 


± 


40 


4910 


± 


53 


5610 it 


83 


Final temperatures (K) 




















Spectroscopy (^12.111 










4900 


± 


100 


5710 it 


100 


0LD, V, and BCv (§13 ■ • 










4970 


± 


154 


5687 it 


130 


Photometry'* 




5207 


± 


108 


4910 


± 


113 


5610 it 


130 


Weighted average 










4917 


± 


67 


5677 it 


68 



^ The number of individual values used to form these weighted averages (from Table fTo| l is indicated in parentheses. The number of 
individual measurements used in each of these weighted averages is indicated in parentheses. All red and near-infrared photometry 
was transformed to the standard system in column 1 as described in the text.^ Color index is slightly out of the range of validity of 
the color/temperature calibration used, but the extrapolation is not considered serious.*^ Uncertainties increased by adding 100 K 



in quadrature to account for possible systematics (see text). 
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TABLE 14 

Estimates of the rotational broadening for the 

COMPONENTS OF CaPELLA. 





Da sin i 


tiB sin i 


Source 


(kms-i) 


(km s^^) 



RJ85 
~10-12 
~10-12 

<30 

35 ± 4 
36 ± (1-2) 

36 ± 3 



36 
36 
35 

32.7 ± 3.3 



36.0 ± 1.5 



^ The measurement presumably refers to the hot star. A 10% uncer- 
tainty has been adopted following Ide Medeiros &: Mavoil lll999fl This 
refers to the total line broadening (see text). 



TABLE 15 
Physical parameters of Capella. 



Parameter Primary Secondary 



Mass (M0) 2.466 ± 0.018 2.443 ± 0.013 

q = Mb/Ma 0.9903 ± 0.0036 

a (10® km) 110.06 ± 0.23 

a (AU) 0.7357 ± 0.0015 

Ttorb (mas) 76.67 ± 0.17 

Distance (pc) 13.042 ± 0.028 

Radius (Rq) 11.87 ± 0.56 8.75 ± 0.32 

logg (cgs) 2.681 ± 0.041 2.942 ± 0.032 

Teff (K) 4920 ± 70 5680 ± 70 

Luminosity (L©)'' 79.5 ± 4.8 72.1 ± 3.6 

BCv (mag) -0.334 ± 0.063 -0.098 ± 0.052 

Mv (mag) 0.315 ± 0.017 0.186 ± 0.016 

t^sini (kms-l)'' 5.0 ± 1.0 35.6 ± 1.5 

t^syncsini (kms-l) 3.92 ± 0.19 2.89 ± 0.11 

Age (Myr)= 540-590 



Computed from V, TTorb, and BCv from IFloweil | |1996D . adopting M®_j = 4.732 
(sec footnote I IH .^ Accounts for the broadening produced by macroturbulence (see 
§[8j.'^ The values given represent the range from all models giving good age agree- 
ment between the primary and secondary from Table [16] Tabic [iTl and Table (Ts] 
for Z fs 0.008. 



Herbig fc Spalding (1955) < 15 

FYanklin (1959) 

Boesgaard (1971) 

Ayres & Linsky (1980) 

aicn ct al. (1985) 9±5 

_Griffin & Griffin (1986) 7 ± (2-3) 

Ffekcl et al. (1986) 5±2 

miisong & Xuefu (1987) 14 ± 3 

^rassmcier & Eckel ( 1990) 3 

_Batten et al. (1991) 

Randich et al. (1994) <7 

de Medeiros fc Mayor (1995 1=" 

_Petrov et al. (1996) . .' 8.5 ± 1.0 

This paper*"" 6.5±1.0 
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TAB LE 16 

Best fit stellar evolution models from IgTrardi et al.I l|200(]fl compared with the 

OBSERVATIONS FOR CAPELLA 





Age 


Mass Temperature 


Radius 


Luminosity 




Star 


(Myr) 


(Mg) (K) 


(Rq) 








Case A: Z = 0.008, «ml = 1-68, aov 


= 0.25, no rotation 




Components constrained to have the same age 


; primary on 


the first ascent 




Primary . . 


537.0 


2.457 5016 


11.99 


80.5 


2.22 






-0.5(T +1.4(T 


+0.2(T 


+0.2(7 




Secondary 


537.0 


2.453 5723 


8.55 


69.3 


1.97 






+0.8(7 +0.6cr 


-0.6(7 


-0.8(7 






Case B 


: Z = 0.008, auL = 1-68, Oov 


= 0.25, no rotation 




Components constrained to have the same age; primary on the descending branch 


Primary . . 


537.0 


2.474 5039 


11.94 


81.2 


3.24 






+0.4O- +1.7<T 


-0.1(7 


+0.4(7 




Secondary 


537.0 


2.453 5723 


8.55 


69.3 


1.97 






+0.8CT +0. Go- 


-0.6(7 


-0.8(7 








Measured values 








Primary . . 




2.466 ± 0.018 4920 ± 70 


11.87 + 0.56 79.5 + 4.8 




Secondary 




2.443 ± 0.013 5680 ± 70 


8.75 + 0.32 72.1 + 3.6 




Note. — 


Differences compared to the model, in the 


sense (model 


minus observation), are 


indicated in 


units of the observational uncertainty a. The measured 


values are given 


at the 


bottom of the tabic for reference. 









TABLE 17 

Best fit Claret stellar evolution models compared with the observations for Capella 



Star 


Age 
(Myr) 


Mass Temperature Radius Luminosity 

(Me) (K) (Rq) (Lq) i2C/"C 


C/N 






Mean age 


Case A: Z = 0.008, q^ml = 1-63, Oov = 0.20, no rotation 
r = 553.2 Myr, At = 1.3 Myr = 0.2%; primary on the descending branch 




Primary . . . . 
Secondary . . 


553.9 
552.6 


2.466 4939 12.16 77.8 19.18 
+0.3(7 +0.5(7 -0.4(7 -2.0(7 

2.443 5736 8.49 69.0 85.01 
+0.8(7 -0.8(7 -0.9(7 


0.64 

+ 1.2(7 

3.27 
-0.2(7 


0.48 
2.02 




Mean 


Case B: Z = 0.004, oml = 1-44, ocov = 0.10, no rotation 
age f = 464.0 Myr, At = 4.5 Myr = 1.0%; primary on the first ascent 






Primary . . . . 
Secondary . . 


461.7 
466.2 


2.466 4917 12.20 76.9 19.13 
-0.1(7 +0.6(7 -0.5(7 -2.0(7 
2.443 5694 8.78 71.6 85.01 

+ 0.2(7 +0.1O- -0.1(7 


1.73 
+19.3(7 

3.27 
-0.2(7 


0.64 
0.06 




Mean 


Case C: Z = 0.002, qml = 1-30, Qov = 0.00, no rotation 
age f = 400.0 Myr, At = 4.0 Myr = 1.0%; primary on the first ascent 






Primary . . . . 
Secondary . . 


398.0 
402.0 


2.466 4936 12.17 77.7 18.88 
+0.2(7 +0.5(7 -0.4(7 -2.0(7 

2.443 5646 9.16 75.4 85.01 
-0.5(7 +1.3(7 +0.9(7 


1.93 
+22.7(7 

3.27 
-0.2(7 


0.48 
2.72 




Case D: Z = 0.008, q?ml = 1-40, aov = 0.20, no rotation 
Mean age t = 588.6 Myr, At = 96.5 Myr = 16.4%; primary in the clump 






Primary . . . . 
Secondary . . 


636.8 
540.3 


2.466 4892 12.43 78.2 19.18 
-0.4(7 +1.0(7 -0.3(7 -2.0(7 
2.443 5674 8.92 72.9 85.01 

-0.1(7 +0.5(7 +0.2(7 


0.64 

+ 1.2(7 

3.27 
-0.2(7 


1.24 
0.33 




Mean age 


Case E: Z = 0.008, oml = 1.63, Qov = 0.20, with rotation 
T = 558.0 Myr, At = 1.1 Myr = 0.2%; primary on the descending branch 




Primary . . . . 
Secondary . . 


558.6 
557.5 


2.466 4938 12.16 77.7 19.24 
+0.3(7 +0.5(7 -0.4(7 -1.9(7 

2.443 5731 8.53 69.4 85.01 
+0.7(7 -0.7(7 -0.8(7 


0.64 

+ 1.2(7 

3.27 
-0.2(7 


0.48 
1.56 






Measured values 






Primary . . . . 
Secondary . . 




2.466 + 0.018 4920 + 70 11.87 + 0.56 79.5 + 4.8 27 + 4 
2.443 + 0.013 5680 + 70 8.75 + 0.32 72.1 + 3.6 


0.57 + 0.06 
3.30 + 0.16 





Note. — Differences compared to the model, in the sense (model minus observation), arc indicated in units of the 
observational uncertainty cr. The measured values arc given at the bottom of the tabic for reference. 



TABLE 18 

Best fit TYCHO stellar evolution models compared with the observations for Capella 



Star 


Age 
(Myr) 


Mass 

(Mq) 


Temperature 
(K) 


Radius Luminosity log Ar(Li) 
(Ro) (Le) {logiV(H) = 12) 12C/13C 


C/N 








Case A: Z = 0.45.Zq; primary jiml = 100; secondary gML = li no rotation 
Mean age f = 593.4 Myr, At = 34.2 Myr = 5.8%; primary in the clump 






Primary .... 
Secondary . . 


610.5 
576.3 


2.466 
2.443 


4935 
+0.2(7 

5736 
+0.8O- 


12.15 78.3 1.87 19.60 
+0.5(7 -0.2(7 +2.8(7 -1.8(7 
8.18 64.7 3.01 82.39 

-1.8(7 -2.1(7 -0.6(7 


0.66 

+ 1.5(7 

3.11 

— 1.2(7 


0.35 
8.03 






Case B: Z = 0.45Zq; primary gML = 1; secondary pML = 1, no rotation 
Mean age f = 537.4 Myr, Ar = 77.8 Myr = 14.5%; primary on the descending branch 






Primary .... 
Secondary . . 


498.5 
576.3 


2.466 
2.443 


4989 
+0.8O- 

5736 
+0.8(7 


11.91 78.7 1.94 18.96 
+0.1(7 -0.2(7 +3.2(7 -2.0(7 
8.18 64.7 3.01 82.39 

-1.8(7 -2.1(7 -0.6(7 


0.69 
+2.0(7 
3.11 

-1.2(7 


1.01 
8.03 










Measured values 






Primary .... 
Secondary . . 




2.466 ±0.018 
2.443 ± 0.013 


4920 ± 70 
5680 ± 70 


11.87 + 0.56 79.5 + 4.8 1.3 + 0.2 27 + 4 
8.75 + 0.32 72.1 + 3.6 3.2 + 0.3 


0.57 + 0.06 
3.30 + 0.16 





Note. — Differences compared to the model, in the sense (model minus observation), are indicated in units of the observational uncertainty 
a. The measured values are given at the bottom of the table for reference. 
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TABLE 19 

Critical times from tidal theory, according to models by 
Claret 

T"sync,A Tsync,B T"circ Consistency with 
Tidal theory (Myr) (Myr) (Myr) observations^ 

Case A: Z = 0.008, ajviL = 1.63, ctov = 0.20, no rotation 
f = 553.2 Myr; primary on the descending branch 

Zahn 545.8 559.8 563.6 yes/yes/no 

Tassoul 540.8 554.6 547.0 yes/yes/yes 

Case B: Z = 0.004, qml = 1.44, Oov = 0.10, no rotation 
f = 464.0 Myr; primary on the first ascent 

Zahn 464.5 471.0 469.9 no/yes/no 

Tassoul 460.3 466.7 466.7 yes/yes/no 

Case C: Z = 0.002, qml = 1.30, aov = 0.00, no rotation 
f = 400.0 Myr; primary on the first ascent 

Zahn 401.8 407.4 406.4 no/yes/no 

Tassoul 396.3 401.8 400.9 yes/yes/no 

Case D: Z = 0.008, «ml = 1-40, ctov = 0.20, no rotation 
T = 588.6 Myr; primary in the clump 

Zahn 537.0 540.8 542.0 yes/no/yes 

Tassoul 532.1 539.5 538.3 yes/no/yes 

Case E: Z = 0.008, qml = 1.63, Oov = 0.20, with rotation 
f = 558.0 Myr; primary on the descending branch 

Zahn 561.0 568.8 574.1 no/yes/no 

Tassoul 555.9 563.6 562.3 yes/yes/no 

Note. — rsync,A and rsync.B are the synchronization times, 
Tcirc is the circularization time, and r represents the mean age of 
the components from Table [TTl 

Consistency with the observations for each of rsync.A, isync.B, 
and Tcirc is 'y^s' if Tsync.A < TA (synchronized primary), rsync.B > 
TB (asynchronous secondary), and rdrc < f (circular orbit). 



TABLE 20 
Coronal abundances for, Capella. 





PIP a 














Aud03 


Aud03 


ArgOS 


Arg03 






Element 


(eV) 


Dra94 


Bri96 


Whi96 


Bau96 


Fav97 


BriOl 


APEC 


MEKAL 


2000 


2001 


Gu06 


Mean 


[C/H].... 


11.3 














-0.39 


-0.36 


-0.09 


-0.12 


-0.20 


-0.23 


[N/H].... 


14.5 


< +0.06 












-0.06 


-0.05 


+0.08 


+0.09 


+0.17 


+0.05 


[U/HJ.... 


13.6 


-0.79 


-0.13*= 


—0.07 








-0.44 


-0.39 


—0.30 


-0.26 


-0.25 


—0.33 


[Ne/H]<^ . . 


21.6 


< -0.41 




-0.16 






-0.58 


-0.22 


-0.42 


-0.30 


-0.30 


-0.15 


-0.30 


[Mg/H] . . . 


7.6 


+0.01 




-0.17 






+0.13 


+0.05 


+0.02 


-0.05 


-0.02 


+0.08 


+0.01 


[Si/H] .... 


8.2 


-0.07 


+0.30 


-0.57 






+0.11 


+0.08 


+0.02 


-0.06 


-0.02 


+0.06 


-0.02 


[S/H] 


10.4 


-0.26 


-0.06 


-0.47 






-0.08 


-0.69: 


-0.94: 






-0.36 


-0.25 


[Ar/H]... 


15.8 




+0.54: 










< -0.06 


< -0.56 






-0.31 


-0.31 


[Ca/H] 


6.1 














-0.54: 


< -1.00 






+0.04 


+0.04 


[Fe/H] . . . . 


7.9 


-0.17 


+0.11 


+0.00 


-0.53: 


+0.00 


+0.03 


-0.04 


-0.05 


+0.01 


+0.01 


+0.07 


+0.00 


[Ni/H] . . . . 


7.6 




+0.26 


-0.17 






-0.42: 


-0.30 


+0.04 


+0.16 


+0.20 


+0.12 


+0.04 


Note. — 


Individual determinations have 


all been rc 


ferred to 


the solar 


abundances 


of IGrevesse & Sauval! 


l|1998|)- 


Upper limits and 


unecrtain 



v alues marked with ':' are not considered in any averag es. Sources of abundance determina tions are as follows: Dra94 ^ iDrake et~an f l994'); Bri96 
^ iBrickhousi U996D: W hi96 ^ 'White? ^1990); Bau96 ^ iBauer Brc gman' (H^gef); Fav97 ^'Fa vata eTal.| (119971) : BriOl ^ iBrickhousec t al. (200j); 
Aud03 — |Audar^^t~a i. (2003), which supersedes the determinations by Audard ct al. (2001). Results listed arc from two independent anal yses 
(APE C, MEKAL); Arg03 ^ Argiroffi ct al.. (^200^. Results listed are from two independent sets of observations (2000, 2001); Gu06 ^ iGu etall 
booel) . Results listed are from a combined RGS+HETGS analysis. 



^ First ionization potential. The separation between low FIP and high FTP is made at 10 cV. Sulfur is near the dividing line, and is considered 
neither low nor high.^ Weighted average of two determinations."^ Strictly speaking, the Nc abundance refers to the coronal value in the Sun, rather 
than the photospheric value, since there are no Ne lines in the solar spectrum. 



